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Sampling  Procedures  and  Tables  for  Life  and  Reliability  TR-3 

Testing  Based  on  the  Weibull  Distribution 
(Mean  Life  Criterion) 


Quality  Control  and  Reliability 


The  content  of  this  technical  report  v&s  prepared  on  behalf 
of  the  Office  of  the  Assistant  Secretary  of  Defense  (installations 
and  Logistics)  by  Professors  Henry  P.  Goode  and  John  H.  K.  Kao  of 
Cornell  University  through  the  cooperation  of  the  Office  of  Haval 
Research.  It  was  developed  to  meet  a  growing  need  for  the  use  of 
mathematically  sound  sampling  plans  for  life  and  reliability  test¬ 
ing  where  the  Weibull  Distribution  adequately  approximates  the  test 
data. 
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FOREWORD 


This  technical  report  presents  a  proposed  acceptance-sampling 
procedure  together  with  tables  of  sampling  plans  for  life  and  relia¬ 
bility  testing  vhen  t ha  underlying  life  distribution  of  items  is  of 
the  Welbull  form.  Cm  study  upon  which  this  report  is  based  was  done 
under  a  Cornell  University  contract  sponsored  by  the  Office  of  Naval 
Research.  Part  of  the  material  was  presented  at  the  Seventh  Rational 
Symposium  on  Reliability  and  Quality  Control  and  is  published  In  the 
Proceedings1  of  this  symposium.  Another  part  of  the  report  was  pre¬ 
sented  at  the  1961  Annual  Convention  of  the  American  Society  for  Quality 
Control  and  was  published  In  the  Transactions2  of  that  convention. 

The  procedures  and  plans  are  for  use  vhen  inspection  of  the  sample 
Items  Is  by  attributes  with  the  life  test  truncated  at  seme  specified 
time.  Lot  quality  Is  evaluated  in  terms  of  mean  Item  life. 

A  set  of  conversion  factors  has  been  provided  frem  which  attribute 
sampling-inspection  plans  of  any  desired  form  may  be  designed  for  the 
Welbull  model  or  from  which  the  operating  characteristics  of  any  given 
plan  may  be  determined.  In  addition,  a  comprehensive  collection  of 
Welbull  sampling-inspection  plans  has  been  designed.  Also  included 
are  tables  of  ratios  for  adapting  the  Military  Standard  105 plans  to 
life  testing  and  reliability  applications.  In  all  three  of  these 
elements  of  the  study  and  the  report,  the  exponential  model  has  been 
included  as  a  special  case  of  the  Welbull  distribution.  Both  the 
procedure  and  plans  are  for  product  for  which  the  value  for  the  Wel¬ 
bull  shape  parameter  is  known  or  can  be  assumed.  Conversion  factors 
and  tables  are  provided  for  a  vide  range  of  values  for  this  parameter. 

These  proposed  techniques  represent  generalizations  of  procedures 
and  plans  developed  for  the  exponential  case  by  Sobel  and  Tischendorf4 
and  by  Epstein**.  Related  work  has  been  done  by  Gupta  and  Groll°  for 
the  gamma  form  of  life -length  distribution.  Recently,  an  interim 
handbook  of  sampling  procedures  and  tables  based  on  the  exponential 
distribution  has  been  published  by  the  Office  of  the  Assistant 
Secretary  of  Defense  (installations  and  Logistics)  based  in  large 
part  on  Professor  Epstein's  work  .  It  is  hoped  that  the  plans  and 
procedures  provided  in  this  report  will  serve  as  a  useful  supplement 
to  this  work. 


Henry  P.  Goode 
John  H.  K.  Kao 

Department  of  Industrial  and  Engineering 
Administration 

Sibley  School  of  Mechanical  Engineering 
Cornell  University,  Ithaca,  New  York 
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SECTION  1 

introduction 

1.1  The  Welbull  Model  In  Life  and  Reliability  Testing. 

Hie  Welbull  distribution  is  a  statistical  nodel  that  has  been  found 
to  be  realistic  for  application  in  life-length  and  reliability  testing  for  many 
mechanical  and  electronic  components.  It  is  a  three -parameter  model  that  takes 
the  form  of  a  general  failure  distribution  in  which  the  failure  rate  need  not 
necessarily  be  constant. 

Further  discussion  of  this  distribution  as  a  life  or  failure  model  will  be 
found  in  Appendix  A  of  this  report.  Reference  may  also  be  made  to  a  paper  by 

O  » 

Kao  .  The  exponential  distribution,  which  is  widely  used  as  a  statistical 
model  in  reliability  analysis  and  inspection,  can  be  considered  as  a  special 
form  of  the  Welbull,  one  which  requires  the  assumption  of  a  constant  failure 
rate. 

The  three  Welbull  parameters  are  (l)  a  "scale"  parameter,  which  is  com¬ 
monly  symbolized  by  the  letter  a,  (2)  a  shape  parameter,  for  which  0  is  the 
symbol  commonly  used,  and  (3)  a  location  parameter,  designated  by  7.  Hie 
general  nature  of  the  Welbull  model  may  be  observed  from  the  density  functions 

plotted  in  Figure  1.  For  use  of  the  procedures  and  plans  covered  by  this 

report,  only  values  for  the  shape  parameter, 0,  and  for  the  location  parameter, 
7,  need  be  determined;  the  value  for  the  "scale"  parameter,  a,  need  not  be 
known. 

From  Figure  1  it  will  be  noted  that  the  form  of  the  distribution  depends 
on  the  shape  parameter,  0.  In  this  figure  the  Welbull  probability  density 
function  has  been  plotted  for  various  values  of  0.  A  plot  for  0  ■  1,  the 

exponential  case,  has  been  included  for  reference.  This  initial  set  of  Welbull 

sampling  plans  is  for  product  for  which  the  value  of  this  parameter  is  known 
or  can  be  assumed  to  approximate  some  given  value.  Conversion  tables  and 
sampling  plans  are  provided  for  nine  values  for  0  ranging  from  l/3  to  5. 

A  relatively  small  number  but  a  broad  range  of  values  for  0  was  selected 
for  this  initial  study.  Ihe  principal  objectives  were  to  develop  practical 
methods  and  techniques  and  to  explore  the  effect  of  differences  in  value  for 
this  parameter  (which  is  the  key  one  for  the  Welbull  distribution).  As  the  use 
of  this  distribution  as  a  statistical  model  Increases,  additional  conversion 
tables  and  sampling  inspection  plans  may  be  constructed  for  intervening  values 
for  0,  particularly  in  the  widely  encountered  region  ranging  from  l/2  to  2. 
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In  May  applications  the  shape  parameter,  p,  nay  be  known  for  the  product 

In  question,  from  past  analysis  of  life  testing  results,  it  may  b<s  established 

that  some  value  of  known  magnitude  may  be  expected  regularly  and  so  may  be  used 

In  sampling  inspection  procedures,  for  example,  for  a  certain  class  of  electron 

q 

tubes  of  receiving  type,  Kao  has  found  from  study  of  approximately  2,000  tubes 
of  a  variety  of  types  and  applications  that  a  value  of  1.7  may  be  appropriate. 
For  ball  bearings,  Lieblein  and  Zelen^found  a  mean  value  of  1.51  with  50#  of 
approximately  5*000  bearings  tested  having  P  in  the  interval  1.17  to  1.7^> 

for  products  for  which  the  value  for  p  is  not  known,  this  parameter  must 
be  estimated  using  failure  data  from  past  inspection  and  research.  A  graphical 
method  for  estimation  will  be  found  in  Appendix  B  of  this  report.  Also,  a 
simple  and  practical  mathematical  procedure  is  currently  being  developed  and 
should  be  available  soonH  further  information  on  estimation  may  also  be  found 
In  papers  by  Kao^^^. 

Also,  as  Indicated,  In  application  of  these  plans  and  procedures,  the  value 
for  7,  the  location  parameter,  must  likewise  be  known  or  assumed  to  approximate 
some  given  value.  For  many  applications  (for  products  for  which  there  Is  risk 
of  failure  immediately  after  manufacture  or  after  being  put  to  use),  a  value  of 
tero  for  7  can  usually  be  assumed.  Methods  for  estimation  of  7  when  estimates 
must  be  made  will  be  found  In  the  appendix  and  In  the  references  cited  above. 

1.2  Form  of  Acceptance  Criteria 

for  the  sampling- inspect  lor.  plans  and  procedures  presented  in  this 
technical  report,  lot  or  product  quail t/  is  evaluated  In  terms  of  mean  item 
life,  p.  Work  is  under  way  on  the  development  of  related  plans  and  procedures 
for  which  the  evaluation  will  be  for  failure  rate  and  for  reliable  life  as 
criteria. 

Inspection  of  sample  Items  is  by  attributes;  all  that  must  be  done  In 
testing  la  to  note  the  ntmber  of  Items  that  have  failed  and  to  do  so  without 
reference  to  the  specific  length  of  life  for  each.  Testing  of  the  sample  Items 
is  stopped  or  truncated  at  the  end  of  seme  preassigned  test  period  length,  t. 

Specifically,  the  following  acceptance-sampling  procedure  for  life  testing 
has  been  assumed: 

1.  Select  at  random  a  sample  of  n  items  from  the  lot. 

2.  Place  the  sample  Items  on  life  test  for  some  preassigned  period  of 
t  time  units. 

3.  Denote  by  y  the  number  of  failures  observed  prior  to  time  t. 
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4.  Accept  the  lot  if  y  i  c,  a  specified  acceptance  number;  if  y  >  c, 
reject  the  lot. 

Curtailed  inspection  for  submitted  lots  prior  to  t  is  possible  for  the  rejection 
of  the  lot  since  it  is  possible  to  observe  (c  +  l)  failures  before  time  t. 

Note  that  this  acceptance  procedure  is  the  same  as  that  specified  for  the 
MIL-STD-105B  sampling  plans  and  other  collections  of  attribute  plans  vith  the 
exception  of  the  introduction  of  a  testing  truncation  time,  t.  It  i6  also  pos¬ 
sible  (as  for  the  105B  plans)  to  employ  double  or  multiple  sampling  instead  of 
single  sampling  as  described  above  and  by  so  doing  reduce  the  average  number  of 
items  at  p'  *  ACJL  that  must  be  put  on  life  test.  However,  the  "economy" 
achieved  is  at  the  expense  of  longer  elapsed  testing  time. 


SECTION  2 

THE  BASIC  CONVERSION  FACTORS 


F 

2.1  The  Function  of  the  Factors. 

The  probability  of  acceptance  for  a  lot,  P(A),  under  plane  of  the 
above  fora  depende  on  the  probability,  p',  of  item  life  being  less  than  (or 
equal  to)  the  teat  truncation  tloe,  t.  For  cases  for  vhich  0  la  known  and 
vith  time,  t,  preaasigned,  p'  la  thus  a  function  of  mean  item  life,  u,  only. 

The  operating  characteristics  of  any  specified  plan  thus  depend  only  on  t  and  pt. 
In  order  to  provide  tables  for  general  use  in  the  design  or  evaluation  of  plana 
for  any  application  rather  than  working  in  terms  of  specific  values  for  t  and 
U,  the  dimensionless  ratio  t/n  is  used.  In  application  of  the  plans  or  tables  to 
any  application,  a  conversion  between  the  ratio  and  specific  t  and  u  values  Is 
extremely  easy  to  make. 

A  set  of  conversion  tables  has  been  computed  to  provide  for  the  Weibull 
distribution  the  connection  between  the  dimensionless  quantity  t/u  and  p' 

(Tables  1  and  2).  With  these  tables,  acceptance- sailing  plans  of  desired  form 
can  be  designed  or  evaluated  using  attribute  sampling  theories  and  practice. 

For  cases  for  which  the  lot  size,  N,  is  large  in  relation  to  the  sample 
size,  n,  the  number  of  failures  prior  to  t  approximates  the  binomial  distribu¬ 
tion  with  parameters  n  and  p',  where  p'  is  defined  as  the  area  under  the  life- 
length  distribution  up  to  t.  The  probability  of  acceptance  P(A)  depends  on  the 
cumulative  number  of  failures  prior  to  time  t.  This  probability  is  given  by 

P(A)  -P(y  <  c)  ■  Z  (?)  p'y  (l-p')n’y  (1) 

y=*o 

The  binomial  distribution  has  been  used  for  the  sampling  plans  given  in  this 
report  except  for  cases  for  which  the  sample  size  is  relatively  large.  For 
these,  the  Poisson  distribution  has  been  used  as  an  approximation  to  the 
binomial.  The  probability  of  acceptance  for  the  Poisson  is  given  by 

P(A)  ■  P  (jt  I  c)  ■  l  (nl>'>y  (2) 

y=o  y! 

An  important  use  for  the  conversion  tables  provided  in  this  report  is  in 
the  adaptation  of  the  MIL-STD-105B  plans  to  reliability  and  life-testing  applica¬ 
tions.  In  describing  the  operating  characteristics  of  these  plans,  the  quality 
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of  submitted  lots  Is  measured  in  terns  of  p',  the  per  cent  defective.  With  the 
conversion  factors  this  form  of  description  may  be  converted  directly  to  measure¬ 
ment  in  terms  of  the  t/p  ratio.  With  this  conversion  the  10f3  plans  may  be 
cataloged  for  appropriate  choice  in  reliability  applications.  Alternatively!  if 
some  105B  plan  has  been  selected!  its  operating  characteristic  curve  may  be 
determined  in  terms  of  the  t/p  ratio,  or  if  the  testing  time,  t,  has  been  speci¬ 
fied!  In  terms  of  the  lot  mean,  4.  An  example  employing  such  a  conversion  is 
shown  later  in  this  report.  It  should  also  be  noted  that  with  the  matching  plans 
provided  in  the  105B  collection!  the  options  of  double- sampling  and  multiple- 
sampling  are  also  available.  Ihe  sample  sites  and  acceptance  numbers  listed  may 
be  used  and  the  established  procedures  for  employing  this  form  of  sampling  in 
attribute  inspection  may  be  followed.  Special  tables  of  conversion  ratios, 
together  with  procedures  for  their  use,  for  directly  applying  the  MIL-STD-105B 
plans  to  reliability  and  life-testing  are  included  in  Section  5  of  this  report. 


{ 
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itation  of  the  Conversion  Tables. 


The  probability  p',  of  an  item  failing  prior  to  acme  testing  time  t 
is  the  value  of  the  cumulative  density  function  at  t.  For  the  Welbull  model 
assumed  in  this  report  [Equation  ( Al6) ],  this  probability  is, 

p'  ■  1  -  exp  (  -(t/iO^I  (3) 

Since  the  mean  of  this  Welbull  distribution  is,  |i  *  i)  T  (j  +  l),  p'  may  be 

0 

re-written  as, 

p'  -  1  -  e::p  {  -  (  £  T  (  i  ♦  1  )  ]P)  (k) 

P  0 


Solving  for  t/p  , 

t/p  ■  [  -ln(l-p* )  J1^  /  r  (  -  ♦  1)  (5) 

0 

Ibis  equation  establishes  the  relationship  between  the  dimensionless  ratio 
t/u  and  p',  the  probability  of  item  life  being  equal  to  or  less  than  t. 

It  may  be  noted  that  for  the  attribute  form  of  sampling  inspection  con¬ 
sidered  here,  only  this  dimensionless  ratio  between  test  time,  t,  and  item 
mean  life,  p,  need  be  of  concern.  Ihe  Welbull  scale  parameter,  r\  has  been 
eliminated.  In  the  mathematics  of  these  plans  and  procedures  it  has  been  assumed 
that  the  Welbull  location  parameter,  7,  has  a  value  of  0  [see  Equation  (Al6)]. 

If  in  application,  however,  7  has  seme  non-zero  value,  all  that  is  necessary  is 
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to  subtract  the  value  for  7  from  the  value  for  t  to  get  tQ,  and  from  the  true 
lot  mean  n,  to  get  uQ«  These  converted  values,  tQ  and  uQ>  are  then  used  for 
all  computations.  Any  solutions  in  terms  of  t  or  uo  can  be  readily  converted 
back  to  real  values  by  simply  adding  the  value  for  7.  This  procedure  for 
handling  the  location  parameter  will  be  Illustrated  later  in  Example  3  and  in 
Example  9.  Of  course,  when  7  -  0,  only  the  parameter  p  (or  b,  which  is  l/p) 
must  be  known. 

To  put  this  relationship  equation  (Equation  (5)1  in  a  form  for  which 
numerical  values  for  relationships  may  be  more  easily  computed,  the  following 
change  is  made: 

t/u  «  (-ln(l-p' )  l*5  /  r(b+l)  ■  exp  (  In  (-ln(l-p' )  )  /  r(b+l) 

»  exp  (  b  In  [-ln(l-p')l)  /r(b+l)  .  (6) 

Values  for  the  expression 

In  (  -In  ( 1  -  p ' )  ]  (7) 

were  obtained  from  a  table  of  the  inverse  of  the  cumulative  probability  function 

Ik 

of  e:rtremes  prepared  by  the  National  Bureau  of  Standards  .  This  table  tabulates 
the  function 

y  «  -In  (-In  9  )  .  (8) 

y 

By  substituting  (l  -  p')  for  <5  the  negative  value  of  Expression  (7)  is  obtained. 

y 

Values  for  e  raised  to  this  power  were  read  from  the  National  Bureau  of  Standards 

tables  of  the  exponential  function1^.  Values  for  the  gamma  function,  r(b  +  l), 

l6 

were  obtained  from  a  table  prepared  by  Dwight  . 

A  table  of  values  for  the  per  cent  truncation,  (t/u)  x  100$,  for  various 
values  of  p'  has  been  prepared.  It  is  presented  as  Table  1.  Values  for  p'  range 
from  .010$  to  80$  with  the  tabulated  values  selected  in  accordance  with  a 
standard  preferred  number  series.  For  convenience  in  both  tabulation  and  use, 
both  the  ratio  t/u  and  p'  are  expressed  as  percentages  rather  than  decimal  ratios. 

For  determining  without  interpolation  the  value  for  p'  when  some  rounded 
value  for  the  (t/u)  x  100  ratio  is  given,  the  relatively  simple  task  of  preparing 
a  table  of  p'  has  been  carried  out. 

The  table  of  values  for  p'  for  various  values  for  (t/u)  x  100  is  presented 
as  Table  2.  Values  for  (t/u)  x  100  range  from  .010  to  100.  Again,  the  values 
used  for  tabulation  form  a  preferred  number  series.  With  this  alternate  table 
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available  together  vith  the  basic  original  one  (Table  l),  e  conversion  nay 
readily  be  made  either  vay--fron  (t/u)  x  100  to  p*  or  frcn  p'  to  (t/p)  x  100. 

Also,  It  will  be  noted  that  the  tvo  supplement  each  other  In  that  B  values 
giving  a  compressed  range  of  figures  in  one  table  give  an  expanded  range  In  the 
other.  This  allows  for  somewhat  more  precise  interpolation  in  conversion.  The 
tvo  together  supply  basic  data  for  the  design  or  evaluation  of  any  life-testing 
and  reliability  sampling  inspection  plan  based  on  the  Welbull  (or  exponential) 
distribution  and  of  an  attribute  form.  For  General  information,  the  relationship 
between  the  (t/u)  x  100  percentage  and  p'  ad  given  by  these  tvo  tables  has  been 
plotted  in  Figure  2  for  each  of  the  various  B  values. 

2.3  Examplen  of  Application. 

Example  (l) 

One  form  of  application  that  should  be  of  considerable  use  is  that  of 
evaluating  the  quality  protection  afforded  by  a  proposed  or  existing  attribute 
acceptance -sampling  plan.  A  possibility  of  iosediate  interest  is  the  use  of  a 
plan  from  the  MIL-STD-105B  Tables. 

Suppose,  for  example,  that  a  105B  plan  vith  an  Acceptable  (feiallty  Level 
(AQL)  of  2.5^  and  Sample  Size  Letter  J  has  been  proposed  for  use.  Reference 
to  Table  IV-A  of  the  105B  Standard  shows  that  for  single  sampling  a  sample  size 
of  75  items  and  an  acceptance  number  of  4  is  specified.  Suppose  life  testing 
time  is  to  be  80  hours  with  simply  a  count  made  of  the  test  items  falling  by  the 
end  of  that  period.  From  Inspection  experience  with  the  product  to  which  the 
plan  la  to  be  applied.  It  seems  most  appropriate  to  assume  a  Welbull  distribution 
with  a  value  for  B  of  1  2/3.  The  lot  size  will  be  relatively  large  compared  to 
the  sample  size  of  75  so  binomial  probabilities  for  sample  items  can  be  assumed. 
Actually,  Table  III  of  MIL-STD-105B  specifies  that  the  lot  size  should  be  from 
1300  to  3200,  501  to  800,  and  101  to  300  for  Inspection  Levels  I,  II,  and  III 
respectively. 

The  first  step  Is  to  determine  the  probability  of  acceptance,  P(A),  for 
various  values  of  p1.  These  probabilities  can  easily  be  obtained  from  any  one  of 
the  readily  available  tables  of  the  emulative  binomial  terms  or  tables  of  the 
Incomplete  beta  distribution.  They  may  also  be  read  from  the  operating  character¬ 
istic  curves  published  as  a  part  of  the  MIL-STD-105B  Tables  (Table  VI).  A  few 
of  these  values  for  this  plan  are  shown  in  the  first  and  second  columns  of  the 
tabulation  below.  Next,  the  first  of  the  conversion  tables,  Table  1,  is  used  to 
obtain  values  for  the  ratio  (t/p)  x  100  for  each  of  the  p*  values,  lhese  table 
values  are  listed  in  the  third  column.  Finally,  using  the  value  for  t  of  80  hours, 
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Figure  2.  The  Relationship  Between  (t./|j)xlOO  and  p1  for  Various  Values  of 


each  of  the  (t/u)  x  100  ratios  are  converted  to  values  for  4.  For  example,  the 
ratio  for  a  p'  of  5^  Is  18.84.  Thus  (80/4)  x  l60  •  18.84  or  u  ■  425  hours.  These 
computations  have  been  made  with  results  as  shown  In  the  last  column  of  the 
tabulation.  One  may  nov  note  that  If  a  lot  Is  submitted  to  this  plan  vhose  mean 
life  is  215  hours,  the  probability  of  its  acceptance  is  only  .01  or  one  in  a 
hundred;  on  the  other  hand,  if  the  mean  life  for  a  lot  is  745  hours  the  proba¬ 
bility  of  acceptance  13  .98.  These  probability  and  mean  life  figures  based  upon 
t  «  80  hours  can  be  plotted,  if  desired,  to  give  the  operating  characteristic 
curve.  (Of  course  similar  OC  curves  for  other  known  values  of  t  may  be  plotted.) 
This  curve  is  the  one  shown  for  3-1  2/3  in  Figure  3* 


Results  of  Computations  -  Example  1 


(in  i) 

P(A) 

(t/u)  x  100 

n 

2 

.98 

10.77 

745 

3 

.92 

13.78 

580 

4 

.82 

16.42 

490 

5 

.68 

18.84 

425 

6  1/2 

.46 

22.15 

360 

8 

.27 

25.20 

315 

10 

.12 

29.01 

275 

12 

.04 

32.58 

245 

15 

.01 

37.63 

215 

To  indicate  the  importance  of  considering  the  shape  of  the  life  density  for 

a  product,  operating  characteristic  curves  for  this  plan  have  been  cjmputed  and 

plotted  in  Figure  3  for  other  selected  values  for  3.  Included  is  a  curve  for  the 

case  in  which  3  equals  1.  This  represents  the  exponential  distribution  widely 

used  as  a  model  in  reliability  and  life-testing  sampling  inspection.  From  these 

curves  it  may  be  noted  that  if  the  underlying  distribution  is  actually  of  the 

Welbull  form  and  the  exponential  is  assumed,  the  actual  operating  characteristics 

of  the  plan  may  differ  very  much  from  those  contemplated.  A  discussion  of  the 

sensitivity  of  statistical  procedures  in  current  use  to  departures  frcm  the 

17 

assumed  exponentiallty  will  be  found  in  a  paper  by  Zelen  and  Dannemlller  . 

It  may  be  noted  in  connection  with  this  example  that  the  MIL-STD-105B  plans 
include  matching  double  and  multiple  sampling  plans.  These  offer  alternative 
possibilities  for  reliability  and  life  testing  applications.  If  incoming  lots  are 
either  quite  good  or  quite  bad  (as  is  commonly  the  case),  substantial  reductions 
in  the  number  of  items  that  must  be  tested  may  be  made.  If  items  are  expensive 
and  if  testing  is  destructive  (as  it  most  likely  vill  be  in  life  testing),  a 
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ure  3.  Operating  Characteristic  Curves  for  Various  Values  of  0 


reduction  In  average  sample  cize  may  be  of  importance.  If  the  teat  period,  t, 
is  relatively  long,  hovever,  the  elapsed  time  required  for  testing  a  second 
sample  (or  subsequent  ones  in  multiple  sampling)  when  such  samples  sure  required 
to  reach  a  decision  may  raise  lfflculties. 

Example  ( 2 ) 

For  a  second  example,  consider  the  case  of  a  manufacturer  who  knows 
that  his  current  production  of  a  certain  component  has  a  mean  life  of  approxi¬ 
mately  52,000  cycles.  Furthermore,  he  has  learned  from  past  experience  with 
life  testing  of  these  components  that  he  can  assume  a  0  value  of  l/2.  A  life 
test  period  of  1000  cycles  seems  Justifiable  and  facilities  are  available  for 
testing  a  sample  of  150  items  from  each  lot.  This  manufacturer  would  like  to 
know  what  acceptance  criteria  to  apply  so  that  virtually  all  lots  will  be  passed 
as  long  as  the  expected  mean  life  of  52,000  cycles  is  maintained.  He  would  aleo 
like  to  know  what  consumer  protection  will  be  afforded.  A  final  question  is 
whether  for  this  application  a  change  to  a  proposed  test  period  of  300  cycles 
and  a  sample  size  of  500  items  would  yield  comparable  or  better  quality  assurance. 

The  first  step  toward  answers  to  these  questions  is  to  compute  the  (t/u)  x 
100  ratio  at  the  mean  life  considered  acceptable.  This  ratio  is  (1000/52,000)  x 
100  or  1.93*  Entering  Table  2  with  this  value  gives  (with  rough  interpolation)  a 
value  for  p'  of  18%  for  a  0  value  of  l/2.  Assuming  a  probability  of  acceptance 
of  .95  is  desired  for  lots  at  the  acceptable  quality  level  of  52,000  cycles  for 
the  lot  mean,  entering  a  table  of  the  cumulative  binomial  distribution  indicates 
that  an  acceptance  number,  c,  of  35  items  gives  this  probability  for  a  sample 
size  of  150  items.  This,  then,  is  the  desired  acceptance  criteria. 

A  simple  measure  of  consumer's  protection  is  to  find  the  lot  mean  value  at 
which  lots  will  likely  be  rejected.  Suppose  a  probability  of  rejection  of  .90 
(of  acceptance  of  .10)  seems  to  be  a  meaningful  figure.  Reference  again  to  a 
binomial  table  indicates  that  for  n  ■  150  and  c  ■  35,  the  probability  of  rejection 
is  .90  at  a  p'  of  approximately  28. Entering  Table  2  with  this  value  gives  a 
(t/n)  x  100  ratio  of  approximately  5-7.  Substituting  a  value  of  1000  cycles  for 
t  in  this  ratio  and  solving  for  u  gives  a  lot  mean  of  17^500  cycles.  This  figure 
for  consumer's  protection  can  be  interpreted  as  follows --since  this  quality 
(u  ■  17,500)  corresponds  to  a  P( A)  ■  .10,  under  this  sampling  plan  (n  ■  150, 
c  ■  35)  on  the  average  90^  of  the  lots  passed  to  the  consumer  will  have  a  mean 
life  of  no  less  than  17,500  cycles.  This  may  or  may  not  represent  adequate 
consumer  protection.  If  it  does  not,  a  plan  with  a  larger  sample  size  must  be 
designed  and  used. 
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An  ansver  to  the  third  question  may  be  found  by  making  similar  computations 
for  an  n  of  500  items  and  a  value  for  t  of  300  cycles.  In  this  case  (t/n)  x  100 
equals  (300/52,000)  x  100  or  .58.  From  Table  2  it  is  found  that  p'  is  approxi¬ 
mately  10$  at  thin  truncation  ratio.  Scanning  a  binomial  table  Indicates  an 
acceptance  number  of  62  vtll  give  a  probability  of  acceptance  of  .95  or  more  when 
the  sample  size  is  500  items.  With  this  sample  size  and  acceptance  number,  the 
probability  of  rejection  is  .90  at  a  p'  value  of  approximately  lU$.  With  this 
value  for  p',  a  (t/n)  x  100  value  of  approximately  115  is  found  from  Table  1. 
Substituting  300  cycles  for  t  in  this  ratio  gives  a  lot  mean  value  of  26,100 
cycles  as  compared  to  17,500  cycles  for  the  first  plan.  Thus  this  combination 
of  sample  size  and  length  of  test  period  gives  better  discrimination  between  good 
and  bad  lots  and  the  consumer  is  therefore  better  protected. 

Example  (3) 

In  this  example,  reference  will  be  made  to  a  case  for  which  the  com¬ 
ponent  life  can  best  be  characterized  by  a  mixture  of  two  Weibull  distributions. 
Kao1^  gives  an  example  of  this  for  the  life  of  electron  tubes.  From  electron 
tube  life  experience,  the  wearout  failures,  i.e.,  drift  of  electrical  properties 
beyond  some  set  limits,  invariably  occur  near  the  latter  part  of  life.  Hence 
the  failures  of  electron  tubes  are  classified  both  as  of  the  wearcut  type  and  as 
of  the  non-wearout  or  catastrophic  type,  each  type  being  represented  by  a  sub¬ 
population  of  the  whole.  In  electronic  terms,  these  failure  types  are  referred 
to  as  electrical  rejects  and  inoperative  rejects  respectively.  The  catastrophic 
(or  inoperative  rejects)  sub-population  is  assumed  to  start  at  time  zero,  i.e., 
the  location  parameter  7^  *  0,  when  the  components  are  first  exposed  to  risks. 

The  wearout  (or  electrical  rejects)  sub -population  is  assumed  not  to  start  until 
some  delayed  period  has  elapsed,  i.e.,  7^  >  0,  since  the  limits  set  on  the  com¬ 
ponent  drift  depend  on  many  factors  such  as  environmental  stress,  maintenance 
policy,  legal  regulations,  and  like  factors.  Since,  in  general,  failures  due  to 
wearout  and  to  non-wearout  reasons  are  identifiable,  it  is  possible  to  treat  the 
two  sub-populations  separately. 

Suppose  that  for  some  application  of  electron  tubes,  the  manufacturer's  past 
experience  indicates  that  the  Weibull  shape  parameter  p^,  associated  with  the 
catastrophic  sub-population  has  a  value  of  l/2,  and  the  shape  parameter  associated 
with  the  wearout  sub-population,  Pg,  has  a  value  of  3  l/ 3 -  Suppose  further  that 
electrical  drift  or  wearout  failure  has  never  been  experienced  prior  to  1000  hours 
of  life.  Under  these  conditions  the  location  parameter  values  will  be: 

7^  =  0  and  “  1000-  Suppose  further  that  the  manufacturer  knows  that  approxi¬ 
mately  2  l/2#  of  the  total  tube  failures  are  of  the  inoperative  type  and  that  the 
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mean  tube  life  for  his  current  production  Is, 

M  =  (.025)  (25,000)  +  (.975)  (11,000)  -  11,325  hrs. 

where  m^  ■  25,000  Is  the  mean  life  of  the  catastrophic  sub -population  and 
Mo  =  11,000  is  that  of  the  vearout  sub-population.  (See  the  appendix  of  the 
paper  by  Kao  J  for  the  derivation  of  this  formula.)  A  life  test  period  of  500 
hours  for  inoperatives  and  of  5000  hours  for  electrical  drifts  are  recommended 
and  acceptance  numbers  -  2  and  Cg  ■  2  for  each  failure  type  are  considered 
satisfactory.  What  ore  the  necessary  sample  sizes  so  that  the  producer's  risk 
is  no  more  than  5%7  Also,  what  would  be  the  consumer's  protection  under  this 
sampling  plan? 

To  answer  these  questions,  the  two  sub -populations  are  treated  separately 
and  are  denoted  by  subscripts  1  and  2  as  done  before  for  the  inoperatives  and 
electrical  drifts  respectively.  For  inoperatives,  (t^/u^)  x  100  ■ 

(500  x  100)  /25,000  *  2.0.  Entering  Table  2  with  this  value  gives  a  value  for 
p'^  of  10.13^  for  a  0  value  of  l/2.  From  a  binomial  table  with  P(A)  k  .95  and 
p'^  *  18.13^  a  value  for  =  5  is  obtained.  Ihe  same  binomial  table  for 
=  5,  *  2  and  P(A)  S  .10,  gives  p'^  ■  75^.  Entering  Table  2  with  this 

value  giveB  (t^/u^)  x  100  =  96.5  and  m1  "  (500  x  100)  /96.5  *  518.24  hours,  a 
value  which  will  be  commented  on  later.  For  electrical  driftB,  7g  must  be  sub¬ 
tracted  from  tg  and  Mg  giving  new  values  for  tg  and  Mg  equal  to  4,000  and  10,000 
respectively.  Hence,  (tg/ug)  x  100  ■  (4,000  x  100)  /l0,000  ■  40.0.  Entering 
Table  2  with  this  value  gives  for  Pg  a  figure  of  3-25^  for  a  0  value  of  3  1/3. 
From  a  binomial  table  with  P(A)  k  .95  and  Pg  ■  1-25%,  it  1b  found  that  ng  ■  25- 
The  same  binomial  table  for  ng  *  25,  Cg  -  2,  P(A)  i  .10  gives  Pg  ■  20%. 

Entering  Table  1  with  this  value  gives  (tg/Ug)  x  100  ®  71.04.  Thus  Mg  * 

(4,000  x  100)  /71.04  -  5,631  hours,  which  upon  reoddlng  7, g  gives  6,631  hours. 
Combining  this  value  of  the  corrected  Mg  with  m^  obtained  for  inoperatives  gives 
the  consumer's  protection  expressed  in  terms  of  a  mean  value  equal  to, 

M  -  (.025)  (518.24)  +  (.975)  (6,631)  -  6,478  hours. 

This  means  under  the  sampling  plan  of  running  a  life  test  for  inoperatives  of 
n  ■  5  and  c  ■  2  for  500  hours  and  another  life  test  for  electrical  drifts  of 
n  ■  25  and  c  ■  2  for  5000  hours,  90%  of  the  lots  passed  to  consimers  will  have 
a  mean  life  of  at  least  6,478  hours.  To  illustrate  the  danger  of  extrapolation 
in  a  mixed  distribution  case,  assume  that  the  second  test  of  5000  hours  duration 
was  net  run  at  all.  Then  the  producer  could  only  base  his  conclusion  upon  the 
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500-hour  test  and  claia  aa  consumers'  protection,  with  90^  confidence,  a  mean 
life  of  at  least  5l8*2U  hours,  a  result  which  iB  altogether  too  modest. 
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SECTION  3 

THE  TABLES  OF  SAMPLING  PLANS 

3-1  Construction  and  Use  of  the  Basic  Tables. 

A  basic  set  of  tables  of  sampling  inspection  plans  has  been  prepared, 
one  table  for  each  of  the  nine  values  for  P  for  which  the  relationship  between 
p'  and  (t/u)  x  100  has  been  established.  These  are  presented  near  the  end  of 
this  report  as  Tables  3a  through  31. 

Each  table  lists  values  for  the  acceptance  number,  c,  and  for  the  minimal 
sample  size,  n,  for  a  variety  of  objective  t /u  ratios.  The  plans  are  designed 
so  that  If  100  times  the  ratio  between  the  test  time,  t,  and  the  mean  life  value 
for  the  lot,  u,  or  (t/u)  x  100  Is  equal  to  or  greater  than  the  selected  column 
value  In  the  table,  the  probability  of  acceptance,  P(A),  will  be  .10  or  less. 
Stated  otherwise,  the  plans  assure  with  90^  confidence  or  more  the  acceptance  of 
lots  for  which  the  (t/u)  x  100  ratio  is  equal  to  or  less  than  the  selected  column 
or  objective  value.  The  ratios  In  the  column  headings  (for  which  the  plans  have 
been  designed)  may  thus  be  considered  in  the  same  way  as  lot  tolerance  per  cent 
defective  (LTPD)  values  are  considered  in  describing  operating  characteristics 
of  the  widely  used  attribute  and  variables  acceptance  plans. 

It  has  been  assumed  that  in  acceptance  Inspection  for  reliability  the 
consumer's  risk  will  be  of  primary  concern.  For  this  reason,  these  plans  have 
been  cataloged  by  P( A)  ■  .10  ratios  vhlch  measure  consumer  protection.  However, 
In  addition  for  each  plan  the  (t/u)  x  100  ratio  is  given  for  which  the  pro¬ 
bability  of  acceptance  is  .95  or  more.  Each  such  P(A)  ■  .95  ratio  value  is  shown 
In  parentheses  under  the  corresponding  srmple  size  number.  These  ratio  values 
may  be  considered  similar  to  acceptable  quality  level  (ktfj)  values  In  Indicating 
the  producer's  risk.  If  the  mean  life  for  the  Items  in  the  lot  Is  such  that  the 
(t/u)  x  100  ratio  is  equal  to  or  lrss  than  the  tabulated  value,  there  Is 
assurance  with  confidence  of  95^  or  more  that  the  lot  will  be  accepted. 

The  two  ratio  values,  one  In  the  column  heading  and  the  other  in  the  body 
of  the  table  In  parentheses  below  the  sample  size  number,  broadly  describe  the 
operating  characteristics  of  each  plan  and  so  form  a  basis  for  making  an  appro¬ 
priate  choice  for  any  acceptance  Inspection  application.  These  values  may  also 
be  used  to  determine  approximately  the  operating  characteristics  of  any 
acceptance  plan  that  has  been  specified  or  that  is  In  use  and  for  vhlch  n  and  c 
match  closely  one  of  the  table  plans.  It  is  easy  to  convert  these  ratios  to 
hours,  cycles  or  some  other  measure  of  lifelength  to  fit  the  product  and  test 
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specif ications  involved.  Thai  will  be  illustrated  by  two  examples  which  follow 
later. 

l8 

In  the  preparation  of  these  plans,  binomial  tables  prepared  by  Grubbs  were 
employed  for  values  for  c  up  to  9  and  for  n  up  to  150.  For  higher  values  of  c  and 
for  values  for  n  up  to  60  or  ao,  the  Pearson  tables  of  the  incomplete  beta-func¬ 
tion  were  used?-^.  Higher  values  of  n  were  determined  by  the  Poisson  approxima- 

20 

tion,  using  a  table  of  np'  values  prepared  by  Cameron  .  The  Poisson  match  was 
checked  and  was  found  close,  even  for  the  smaller  sample  sizes  and  the  larger 
values  for  p'.  The  slight  differences  that  may  exist  in  some  cases  is  on  the 
conservative  side;  the  value  for  n  is  slightly  larger  than  that  theoretically 
require!.  As  the  tables  were  prepared  as  a  part  of  an  exploratory  study,  plans 
showing  extremely  large  sample  sizes  have  been  included  to  indicate  the  order  of 
magnitude  Involved  and  not  with  the  expectation  that  samples  of  this  size  would 
ordinarily  be  used  in  acceptance-inspection  practice. 

3*2  Examples  of  Application  of  the  Tables. 

Example  (4) 

Suppose  an  acceptance  inspection  plan  is  required  which  will  aecure 
with  90^  confidence  a  mean  life  for  items  of  4000  hours  or  more  for  each  lot 
accepted.  Also,  suppose  it  will  be  desirable  to  assure  the  producer  that  if  the 
mean  life  for  items  in  a  lot  is  25,000  hours  or  more,  there  will  be  a  high 
probability  (.95)  that  the  lot  will  be  accepted.  A  test  period  of  400  hour3  for 
the  inspection  of  sample  items  has  been  specified.  Through  past  experience  it 
has  been  determined  that  the  distribution  of  item  life  is  of  the  Weibull  form 
with  the  shape  parameter,  0,  equal  to  approximately  l/2.  Also,  a  value  for 
7,  the  location  parameter,  of  0  can  be  assumed. 

For  these  sampling  plan  specifications,  100  times  the  ratio  of  test  time,  t, 
to  the  mean  life,  ^  1Q,  for  which  a  probability  of  acceptance  of  .10  or  less  is 
desired  is  (400/4,000)  x  100  or  10.  At  the  .95  probability  value  the  ratio  is 
(400/25,000)  x  100  or  1.6.  A  plan  approximately  meeting  these  requirements  may 
be  found  in  Table  jb  which  gives  plans  for  distributions  for  which  P  *  l/2.  The 
column  for  which  (t/u)  x  100  ■  10  is  entered  and  then  scanned  for  the  ratio  value 
1.6  (or  one  close  to  this  value)  among  the  values  in  the  column  listed  In 
parentheses.  This  value  is  found  well  down  in  the  column.  The  corresponding 
sampling  Inspection  plan  specifies  a  sample  size,  n,  of  43  and  an  acceptance 
number,  c,  of  11. 


! 


Example  (5) 

*.  Suppose  that  a  sampling  Inspection  plan  specifies  that  a  random  sample 

of  3000  items  is  to  be  drawn  from  the  lot  and  tested  for  a  period  of  2,480  hours 
If  no  more  than  7  items  fail  before  the  end  of  the  test  period,  the  lot  is  to  be 
accepted;  if  more  than  7  items  do  not  li'e  through  the  test  period,  the  lot  is 
to  be  rejected.  Life  measurements  for  past  inspection  and  research  for  the 
product  to  which  the  plan  is  to  be  applied  indicate  the  distribution  is  of  the 
Weibull  form  with  0  equal  to  approximately  1  2/3.  Also,  a  value  for  7 ,  the 
location  parameter,  of  tero  is  indicated.  The  prospective  user  of  this  plan 
would  like  to  know  what  quality  protection  will  be  given.  Inspection  of  Table 
3d  which  lists  plans  for  0  =  1  2/3  discloses  a  plan  matching  reasonably  well  the 
one  specified,  the  plan  for  which  c,  the  acceptance  number,  is  7  and  n,  the 
sample  size,  is  3,019.  For  this  table  plan  the  (t/u)  x  100  ratio  at  P(A)  »  .10 
is  4.  Substitution  of  the  specified  test  period  length  of  480  hours  for  t  gives 
(UdO/u)  x  100  «  4.  Solving  for  4  gives  12,000  hours  as  the  mean  value  for  item 
life  for  the  lot  for  which  the  probability  of  acceptance  in  .10  or  less.  A 
similar  substitution  for  t  using  the  ratio  at  which  P(A)  ■  .95  gives 
(480/n)  x  100  =  2.1.  Solving  for  m  again  gives  23,000  hours  as  a  lot  mean  value 
for  which  the  probability  of  acceptance  is  .95  or  better.  Ihe  values  for  the 
lot  mean  at  these  two  probability  values  broadly,  but  very  practically,  describe 
^  the  operating  characteristics  of  the  specified  plan. 

3.  3  Some  Points  of  Practice. 

In  the  use  of  these  tables  of  plans,  several  points  of  practice  should 
be  noted.  First,  in  using  the  p1  values  associated  with  values  for  ( t/fi )  x  100 
for  the  Weibull  distribution  to  find  values  for  c  and  n,  the  binomial  probability 
distribution  has  been  used.  Hence  the  size  of  the  lot  should  be  relatively  large 
compared  to  the  size  of  the  sample  for  the  stated  probability  values  to  pre¬ 
cisely  apply.  Second,  if  a  plan  is  net  available  for  which  a  (t/n)  x  100  ratio 
in  the  column  headings  matches  closely  the  desired  ratio,  to  be  conservative,  a 
plan  should  be  chosen  from  the  column  with  the  next  smaller  ratio  heading.  Ihis 
will  assure  w?  *  h  confidence  greater  than  90^  the  specified  mean  life  for 
acceptance.  T£  the  acceptable  quality  level  (the  ratio  or  mean  life  for  which 
P(A)  ■  .95)  must  also  be  guaranteed  and  a  matching  ratio  value  is  not  found  in 
the  selected  column  of  plans,  a  plan  vlth  the  next  greater  value  should  be 
selected.  Lots  equal  to  or  better  than  the  specified  "acceptable  quality"  will 
have  an  assurance  of  greater  than  95^  of  being  accepted.  With  proper  care,  seme 
rough  interpolation  may  be  employed  between  listed  sample  sizes  (either  down  or 
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across  the  table  or  both)  to  find  a  new  plan  having  more  nearly  the  desired 
characteristics.  Finally,  If  a  plan  Is  found  for  vhlch  the  desired  or  given 
ratios  closely  match  but  for  practical  reasons  it  seems  desirable  to  round  off 
the  sample  size  to  the  nearest  number  ending  in  zero  or  five,  such  rounding  off 
should  be  done  to  the  next  larger  size.  This  will  assure  retention  of  the  proba¬ 
bility  values  of  .10  or  less  for  the  ratios  given  in  the  column  headings. 


SECTION  h 

PLAN  DESIGN  Hi  TERMS  OF  QUALITY  LEVEL  RATIOS 

*♦.]  The  Rela  Ion  ship  Between  Acceptable  and  Unacceptable  Lot  Quality. 

Sampling  plans  are  moat  conveniently  cataloged,  selected,  or  designed 

Ln  terms  of  a  producer's  risk  and  a  consumer's  risk.  Some  lot  quality  figure 

will  be  specified  as  satisfactory  and  for  lots  of  this  quality  or  better  the 

probability  of  acceptance  should  be  high,  conventionally  .95  or  more  (the  pro* 

ducer's  risk  or  rejection  small,  .0$  or  leas).  For  the  plans  for  life  testing, 

described  ln  Section  3,  this  specification  will  be  a  lot  mean  life,  4  ^5 »  at 

which  P( A)  2  .95-  Likewise,  an  unsatisfactory  quality  level  will  be  specified 

for  which  the  probability  of  acceptance  will  be  low,  conventionally  .10.  This 

specification  will  be  a  lot  mean  life,  4  ^  at  which  P(A)  i  .10.  (if  other 

values  of  producer's  and  consumer's  risks  must  be  determined  for  these  plana, 

21 

reference  may  be  made  to  a  method  in  a  paper  by  Kao  .  ) 

In  plan  selection  or  design,  one  objective  is  to  find  a  combination  of 

sample  size  and  acceptance  number  which  simultaneously  yields  the  desired  values 

for  both  the  consumer  risk  and  the  producer  risk.  If  working  from  tables  of 

plans,  the  values  for  lot  quality  at  the  two  risk  figures  may  be  found  listed 

ln  the  tables.  In  the  design  of  a  plan,  one  may  cut  and  try  until  a  suitable 

plan  is  found  in  a  manner  suggested  ln  Example  2.  Also,  factors  are  available 

which,  ln  conjunction  with  the  conversion  tables  supplied  here  as  Table  1  and 

20 

Table  2,  enable  a  direct  determination  to  be  made 

A  simple  alternative  solution  for  the  form  of  acceptance  inspection  dis¬ 
cussed  here  is  to  make  use  of  one  of  its  properties,  namely  that  for  a  given 
acceptance  number,  c,  (and  for  a  given  value  for  0)  the  ratio  betwuen  the  lot 
means  at  the  two  risk  values  is  approximately  constant  for  all  values  of  sample 
size,  n.  These  ratios  (or  multipliers)  have  been  determined  as  part  of  the  study 
underlying  this  report  for  values  for  c  ranging  from  0  to  15  for  each  of  the 
various  values  for  0.  They  are  presented  in  Table  4.  The  table  values  are  in 
the  form  of  multipliers  for  finding  4  ^ ,  given  4  or  by  using  the  reciprocal 
of  the  multiplier,  for  finding  4  given  4  ^5 •  That  M  55  (for  which 

P(A)  =  .95)  is  equal  (approximately)  to  4  ^  (for  which  P(A)  ■  .10)  times  the 
appropriate  table  multiplier.  These  multipliers  may  be  used  both  to  assist  ln 
evaluating  the  operating  characteristics  of  seme  given  plan  and  to  assist  in  the 
design  of  a  plan  to  meet  seme  acceptance-inspection  requirement. 
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**.2  An  Example  of  Use  of  the  Table  of  Multipliers. 


Example  (6) 

For  a  certain  purchased  component  the  lot  mean  life  should  be  at  leari 
4,000  hours;  this  value  is  accordingly  chosen  for  u  .  Alec,  the  producer  has 
been  informed  that  lots  whose  mean  life  is  10,000  hours  or  more  are  reasonably 
sure  of  acceptance  through  the  sampling  procedure.  Accordingly,  this  value  is 
to  be  used  for  u  Qt..  A  value  for  0  of  1  can  be  assumed.  A  testing  period,  t, 
of  200  hours  has  been  specified.  Values  for  sample  size,  n,  and  acceptance 
number,  c,  must  be  found  to  meet  these  requirements. 

The  ratio  between  the  two  lot  means,  u  95/^  10'  ie  10, 000/4,000  or  2.5- 
Examination  of  the  table  of  mean  life  multipliers,  Table  4,  under  the  column 
for  0  *  1  indicates  that  an  acceptance  number,  c,  of  10  items  will  give  this 

ratio.  The  (t/u)  x  100  ratio  at  m  ^  is  (200/4,000)  x  100  or  5-  Entering 

Table  2,  the  table  of  p',  with  this  truncation  ratio  value  of  5,  giveB  a  p’  of 
4.-88$.  Reference  to  a  table  of  the  cumulative  binomial  distribution  or  use  of 
the  Poisson  approximation  for  c  ■  10  and  p'  »  .0488  at  P(A)  =  .10  shows  that  a 
sample  size,  n,  of  315  items  meets  the  requirements.  A  check  for  this  solution 

can  be  made,  if  desired.  For  n  *  315>  c  *  10,  and  P(A)  =  .95  the  Poisson 

approximation  indicates  a  p'  of  1.9<S$.  Entering  Table  1,  the  table  for  per  cent 
truncation,  with  this  value  for  p',  a  value  for  (t/u)  x  100  of  approximately 
2.0  is 
is  the 


found.  Solving  for  u  yields  (200/u)  x  100  «  2.0  or  p  Qt.  =  10,000  which 
desired  value. 
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SECTION  5 

AN  ADAPTATION  OF  THE  MIL- STD- 10 5B  PLANS 

5-1  Use  of  the  105B  Plans  for  Life  and  Reliability  Testing. 

To  permit  use  to  be  made  of  the  familiar  Military  Standard  105B  plans 
for  life  and  reliability  tebtlng  application,  the  conversion  factors  described 
in  Section  2  of  this  report  were  employed  to  find  (t/n)  x  100  ratios  for  all 
of  the  plans  in  the  105B  collection.  As  for  the  basic  plans  described  in 
Section  3,  separate  tables  have  been  prepared  for  each  of  a  number  of  selected 
values  for  0,  the  Weibull  shape  parameter.  The  special  case  for  0-1,  which 
Is  the  exponential  case,  has  been  Included. 

The  acceptance  procedure  will  be  the  same  as  that  employed  for  the  basic 
plans  and  as  outlined  In  Section  1.2:  (l)  a  random  sample  of  n  Items  1b  selected 

from  the  lot,  (2)  the  sample  Items  are  tested  for  life  over  some  preassigned 
test  period  length,  t,  (3)  the  number  of  test  items  failing  prior  to  time  t  is 
observed,  (4)  If  the  number  of  test  items  failing  Is  no  more,  than  some  specified 
acceptance  number,  c,  the  lot  Is  accepted  -  if  more,  the  lot  is  rejected. 

Lot  quality  is  evaluated  In  terms  of  mean  item  life,  4.  Both  t  and  11  are 
measured  from  some  reference  time  or  the  value  for  the  Weibull  location  para¬ 
meter,  y. 

The  sample  sixes  and  acceptance  numbers  used  will  be  those  specified  by 
the  MIL-STD-105B  tables.  The  procedure  as  outlined  here  is  for  single  sampling; 
through  simple  and  appropriate  modification,  the  105B  double -sampling  and 
multi  pie -sampling  plans  may  be  likewise  employed.  Example  7  In  a  following 
subsection  of  this  report  discusses  this  possibility.  It  may  be  noted  that  the 
acceptance  procedure  (for  any  form  of  sampling)  is  the  sjkae  as  that  specified 
for  the  MIL-STD-105B  plans  with  the  single  exception  of  the  use  of  a  test 
truncation  time,  t,  (and  the  application,  of  course,  to  life-testing  data). 

As  discussed  In  Section  2,  the  probability  of  acceptance  for  a  lot  under 
the  acceptance  procedure  outlined  depends  solely  upon  the  probability,  p',  that 
an  item  will  fall  before  the  end  of  the  test  period,  t.  If  the  test  truncation 
time,  t,  is  preassigned  and  If  the  value  for  0,  the  shape  parameter,  is  known 
the  probability,  p',  of  failure  is  a  function  only  of  mean  item  life,  This 
fact  makes  it  possible  to  use  this  attribute  acceptance  procedure  to  evaluate 
lots  in  terms  of  mean  item  life;  the  operating  characteristics  for  any  speci¬ 
fied  plan  (in  terms  of  c  and  n)  depend  only  on  t  and  |i. 


o  •  m  -  4 
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It  should  be  noted  again  that  it  is  necessary  to  assume  a  value  for  P, 
the  Welbull  shape  parameter.  For  many  applications  this  value  may  be  known. 

Its  magnitude  may  have  been  determined  for  the  product  in  question  from  past 
life-length  research  data,  from  the  results  of  past  inspection  data,  or  from  some 
other  source.  If  the  value  for  0  is  not  known,  procedures  are  available,  ns 
outlined  in  Appendix  B,  for  estimating  this  parameter  (and  the  7  location 
parameter  also  if  this,  too,  is  necessary). 

5.2  The  105B  Tables  of  Ratios. 

Tables  of  ratios  for  adapting  the  105B  plans  to  reliability  and  life¬ 
testing  application  have  been  prepared  for  each  of  seven  typical  values  for 
0.  These  seven  values  range  from  l/2  to  3-1/3;  covering  the  span  commonly 
encountered  with  industrial  products.  A  table  for  0  ■  1,  which  is  the  exponential 
case,  is  also  included,  lhese  tables  will  be  found  at  the  end  of  this  report 
as  Tables  5a  through  5g. 

In  order  to  allow  for  any  desired  test-time  truncation  value,  t,  and  to 
make  the  plans  available  for  general  use,  the  tables  have  been  prepared  in 
terms  of  the  dimensionless  ratio,  t/u.  Actually,  to  give  more  conveniently 
usable  figures  and  to  work  in  terms  adopted  for  the  105B  plans,  the  ratio  is 
given  in  terms  of  percent;  (t/u)  x  100  is  used.  Each  of  the  105B  plans  is 
cataloged  and  described  in  terms  of  the  (t/u)  x  100  ratio.  These  ratio  values 
are  used  in  the  same  way  as  the  percent  defective  values  are  used  in  the 
selection  and  application  of  the  105B  plans  for  ordinary  attribute  inspection. 

When  applying  the  plans  to  a  specific  life-testing  application,  use  of  the  ratio 
to  convert  from  test  time  to  lot  mean  (or  vice  versa)  will  be  found  quite  easy. 
Examples  of  application  are  giver,  in  Section  5- 3- 

Each  table  lists  for  each  105B  Acceptable  Quality  Level  (AQL)  value  the 
corresponding  (t/u)  x  100  ratio.  These  matched  ratios  will  be  found  in  the 
column  headings  under  each  of  the  respective  105B  Acceptable  Quality  Level 
values  (which  are  in  termB  of  100  p'  the  acceptable  percent  defective)  Each 
ratio  value  gives  for  all  105B  plans  of  the  corresponding  AQL,  a  measure  of 
lot  quality  for  which  the  producer's  risk  or  probability  of  rejection  will  be 
low.  This  risk  of  rejection  will  be  the  same  as  that  encountered  in  the  normal 
use  of  the  105B  plans  for  attribute  inspection.  It  will  be  recalled  that  this 
risk  is  not  a  constant  vsdue  of,  say  .05,  as  in  most  previous  tables  of  accept¬ 
ance  inspection  plans,  but  ranges  from  as  low  as  0.01  to  as  high  as  0.20.  The 
risk  varies  with  the  size  of  the  sample,  which  in  turn  varies  with  the  size  of 
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the  lot.  For  large  lot  sizes  (and  thus  large  sample  sizes)  the  risk  is 
relatively  email;  for  small  lot  sizes  it  is  relatively  large.  The  specific 
risk  value  for  any  plan  of  interest  may  he  obtained  from  the  corresponding 
operating  characteristic  curve  which  will  be  found  included  with  the  105B  tables. 

The  interpretation  of  these  matched  Acceptable  Quality  Level  ratios  for 
life-testing  and  reliability  use  may  be  demonstrated  by  means  of  a  specific 
case.  Assume,  for  example,  that  0  =  1-2/3  and  that  a  105B  plan  with  an  AQL  of 
b.0%  is  to  be  used.  From  the  table  of  ratios  for  0  *  1-2/3,  which  is  Table  5e, 

it  will  be  found  that  the  corresponding  t/n  ratio  value  at  the  AQL  la  16.42. 

Thus  lots  for  which  ( t/p. )  ::  100  =  l6.k2  are  "acceptable"  and  the  probability  of 
acceptance  will  be  high  (the  probability  of  rejection  low).  If  the  test  period, 
t,  is  ,  say,  1000  hours,  (lOO/p)  x  100  =  16.42  or  p  «  6,090  hours;  the  swan  life 

for  the  items  in  the  lot  must  be  6,090  hours  for  it  to  meet  the  105B  acceptable 

quality  level  standards. 

In  the  body  of  each  table  of  conversion  ratios  will  be  found  for  each  105B 
plan  ratios  for  which  the  probability  of  acceptance  is  .10.  These  correspond 
to  lot  tolerance  per  cent  defective  (LTPD)  figures  which  furnish  a  useful 
measure  of  consumer's  protection.  They  represent  unsatisfactory  lot  quality 
values  and  ones  for  which  the  probability  of  accept  uice  is  low.  Unlike  the 
risks  associated  with  the  AQL  which  vary  with  sample  size,  the  risk  at  the  LTPD 
quality  used  for  the  tables  presented  here  is  at  .10  for  all  plans. 

For  the  example  cited  above  in  which  the  AQL  is  4.0 1>,  suppose  the  Sample 
Size  Code  Letter  to  be  used  is  L.  Reference  to  the  table  of  conversion  ratios 
for  0  =  1-2/3  shows  the  LTPD  ratio  to  be  31.  With  the  test  period,  t,  designated 
as  1000  hours,  the  value  for  u,  the  lot  mean  life,  may  be  readily  determined. 
Substitution  gives  (1000/m)  x  100  =  31  or  n  ■  3>220  hours.  Thus  lots  whose  mean 
life  is  3220  hours  or  less  have  a  probability  of  at  most  .10  of  acceptance. 

With  the  use  of  these  complete  tables  of  conversion  ratios  (tables  5®  through 
5g),  suitable  105B  plans  may  be  selected  in  terms  of  either  an  AQL  or  the  LTPD 
(or  both,  if  desired).  If,  instead,  some  .105B  plan  has  been  specified,  its 
operating  characteristics  can  be  evaluated.  Examples  of  such  use  will  be  out¬ 
lined  in  the  material  that  follows. 

As  a  supplement  to  these  tables,  Table  6  has  been  prepared.  This  table 
gives  the  (t/u)  x  100  ratio  at  the  Acceptable  i^iality  Level  for  an  additional 
number  of  values  for  0,  the  Weibull  shape  parameter.  Ratios  at  the  AQL  for  the 
0  values  used  in  Table  5  are  also  included  for  convenience.  As  the  Acceptable 
Quality  Level  supplies  the  operating  characteristic  of  most  interest  in  the 
application  of  105B  plans,  the  ratio  values  in  this  table  may  be  all  that  are 
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necessary  for  many  applications.  Ratios  are  given  for  Seven  additional  values  for  0 
for  which  conversion  factors  were  available.  Also  included  is  a  table,  Table  7, 
of  the  sample  sizes  and  acceptance  criteria  for  the  Military  Standard  105B  plans 
for  single  sampling. 

5-3  Examples  of  Application. 

Example  (7) 

For  a  simple  example  0*  application,  consider  a  receiving  inspection 
case  for  which  incoming  lots  of  a  product  are  tc  be  tested  for  lifelength  by 
sampling.  From  past  experience  with  the  product  it  has  been  determined  that  the 
life  distribution  can  be  expected  to  follow  the  Weibull  form  with  a  value  for  0, 
the  shape  parameter,  of  approximately  1-1/3-  The  value  for  7,  the  location 

parameter,  is  expected  to  be  0.  The  MIL-STD-105B  plans  are  to  be  employed.  A 

test  period  for  the  sample  items  of  200  hours  and  an  Acceptable  Quality  Level 
in  terms  of  percent  defective  (as  used  in  the  Standard)  of  1.5#  have  been  more 
or  less  arbitrarily  selected  for  use.  The  size  of  incoming  lots  is  5,000  items. 
Inspection  Level  II,  the  one  for  normal  use,  will  seemingly  be  appropriate. 

Single  sampling  is  to  be  employed.  The  acceptance  procedure  for  the  above  condi¬ 
tions  and  the  resulting  operating  characteristics  must  be  determined. 

Reference  to  Table  III  of  MIL-STD-105B  shows  that  for  a  lot  size  of  5,000 
items  and  for  Inspection  Level  II,  Sample  Size  Code  Letter  M  is  designated  for 
ordinary  inspection.  Reference  is  next  to  Table  IV-A  of  the  Standard,  the  master 
table  for  normal  6 ingle -sampling  inspection  or  to  Table  7  of  this  report.  Here 

it  will  be  found  that  for  Sample  Size  Code  Letter  M  and  for  an  Acceptable  Quality 

Level  of  1.59&,  the  sample  size  is  225  items,  the  acceptance  number  is  8  Items, 
and  the  rejection  number  is  9-  The  acceptance-rejection  procedure  will  thus  be 
the  following:  (a)  draw  at  random  from  the  submitted  lot  a  sample  of  225  items 
and  place  them  on  life  test  for  200  hours,  (b)  determine  the  number  of  items  that 
have  failed  by  the  end  of  this  test  period,  (c)  if  the  number  failing  is  8  or 
less,  accept  the  lot.  If  the  number  is  9  or  more,  reject  it. 

The  operating  characteristics  of  this  plan  can  be  determined  from  information 
included  in  the  tables  of  (t/p)  x  100  ratios  included  as  part  of  this  report. 

For  this  example  reference  will  be  to  Table  5<1,  the  table  of  ratios  for  0  »  1-1/3. 
Examination  of  the  two  lines  of  Acceptable  Quality  Level  values  across  the  top 
of  this  table  shows  that  for  an  Acceptable  Quality  Level  in  terms  of  p'  {$)  of 
1.5,  the  Acceptable  Qiallty  Level  in  terms  of  (t/p)  x  100  is  4.69.  With  this 
latter  ratio,  and  with  the  value  for  t,  the  test  period  length  of  200  hours,  the 
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value  for  u,  the  mean  Item  life  for  the  lot,  can  be  determined.  Thus: 

(t/n)  x  100  -  4.69  (aql) 

(200/m)  ::  100  *  4.69 
M  -  4,260  hours. 

One  now  knows  that  the  operation  of  the  plan  Is  such  that  If  the  mean  Item  life 
for  the  lot  la  4,2o0  hours  or  more  the  probability  that  It  will  be  accepted  la 
high.  (A  rough  value  for  this  probability  may  be  found  from  the  operating 
characteristic  curves  In  the  Military  Standard.  For  an  AQL  of  1.5  and  for  Code 
Letter  M  one  may  note  the  probability  is  approximately  .97. )  Ihe  Acceptable 
Quality  Level  is  thus  4,260  hours. 

Ihe  ability  of  the  plan  to  protect  the  consumer  may  be  measured  by  the  lot 
mean  life  for  which  the  probability  of  acceptance  is  low.  Ihe  tables  of  (t/u)  x 
100  ratios  Included  In  this  report  Include  ratios  at  the  Lot  Tolerance  Percent 
Defective  (LTPD)  quality  level,  the  level  at  which  the  prbbability  of  acceptance 
is  .10.  Ihese  ratios  will  be  found  In  the  body  of  the  tables.  Reference  to  the 
same  table,  Table  5<1  for  6  ■  1-1/3*  gives  an  L1PD  ratio  of  13  for  Sample  Size 
Code  Letter  M  and  an  AQL  of  1.5 •  Ccxiputatlons  similar  to  those  previously  made 
give: 

(t/u)  x  100  ■  13  (LTPD) 

(200/m)  x  100  -  13 
M  -  1,540  hours. 

One  now  knows  that  If  the  mean  life  for  Items  in  a  submitted  lot  Is  1,540  hours 
or  less,  the  probability  of  It  being  accepted  Is  at  most  .10;  the  probability  of 
its  rejection  is  at  least  .90. 

Ihe  operating  characteristics  In  hours  as  computed  above  apply  also  (with 
the  same  values)  If  a  double -sampling  or  a  multiple -sampling  plan  for  the  same 
Sample  Size  Code  Letter  and  AQL  value  is  employed  Instead  of  a  single-sampling 
plan.  For  double -sampling  In  this  application,  the  data  for  the  plan  will  be 
found  in  Table  IV-B  of  the  Military  Standard.  Ihe  first  soniple  size  vlll  be  150 
items.  Ihese  sample  items  would  be  tested  for  200  hours.  If  5  or  fewer  Items 
failed  within  this  time,  the  lot  would  be  accepted;  If  14  or  more  failed,  It  would 
be  rejected.  If  from  6  to  13  failed,  a  second  sample  of  300  Items  would  be 
selected  and  tested  for  200  hours.  If  the  total  number  failing  (in  the  first  and 
second  samples  combined)  is  13  or  less,  the  lot  would  be  accepted,  if  it  Is  14  or 
more,  the  lot  would  be  rejected. 

Example  (8) 

For  a  second  example,  consider  an  acceptance-inspection  by  sampling 
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application  for  which  the  following  achievement*  are  desired:  (a)  If  the  mean 
item  life  for  the  lot  is  20,000  hours  or  more  the  probability  of  acceptance  will 
be  high.  Lots  of  this  mean  life  or  greater  are  considered  "acceptable."  (b)  If 
the  mean  item  life  for  the  lot  is  6,000  hours  or  less,  the  probability  of 
acceptance  will  be  low,  that  is  .10.  A  test  period  of  $00  hours  will  be  employed. 
It  is  expected  that  the  item  life  distribution  will  be  of  the  Weibull  form  with  a 
value  for  3,  the  shape  parameter  of  3A'  She  value  for  the  location  parameter, 

7,  is  zero. 

Hie  (t/u)  x  100  ratio  at  the  A<&  will  be 
(500/20,000)  x  100  or  2.5  . 

The  (t/u)  x  100  ratio  at  the  L1TD  will  be 
(500/6,000)  X  100  or  8.3  . 

With  these  values,  Table  5b  giving  the  (t/u)  x  100  ratios  for  3  *  3A  nay  be 
scanned  to  determine  the  appropriate  MIL-STD-105B  plan.  One  may  note  in  this 
table  that  an  AQL  of  6.5  (in  percent  defective)  corresponds  to  a  (t/u)  x  100  ratio 
of  2.29.  lhis  is  the  closest  value  available  for  the  desired  ratio  value  of 
2.5.  Next,  the  column  under  the  6.5  AQL  value  heading  may  be  scanned  to  find  a 
close  approximation  tc  the  desired  value  of  8.3  for  the  LTPD.  A  value  of  8.2, 
which  is  reasonably  close,  is  found  corresponding  to  Sample  Size  Letter  K. 

Thus  any  MIL-STD-105B  plan  with  Sample  Size  Code  Letter  K  and  with  an 
Acceptable  Quality  Level  of  6.5  will  give  approximately  the  desired  operating 
characteristics  for  the  specified  test  period  of  500  hours.  For  single  sampling, 
for  example,  the  sample  size  will  be  110  items  and  the  acceptance  number  12  as 
indicated  in  the  MIL-STO-105B  tables. 

Example  (9). 

The  procedure  to  be  followed  for  esses  in  which  the  Weibull  location 

parameter,  7,  is  not  zero  but  is  of  some  other  known  value  may  be  illustrated  by 

outlining  a  third  example.  The  method  to  be  followed  in  this  case  is  to  simply 

subtract  the  value  for  7  from  the  value  for  t,  the  test  time,  to  get  tQ,  and 

from  u  to  get  u  ■  Biese  transformed  values  t  and  Uj  are  then  used  for  all 
o  00 

(t/u)  x  100  computations.  The  solution  obtained  in  terms  of  t  and  uQ  can  then 
be  converted  back  to  original  values  by  simply  adding  the  value  for  7  to  each. 

Consider,  for  example,  an  application  for  which  a  single- sampling  plan  with 
n  equal  to  35  and.  c  equal  to  1  has  been  specified.  Ibis  corresponds  to  a  plan 
with  Sample  Size  Letter  H  and  an  AGE,  (in  terms  of  p’)  of  1.5$  in  the  105B 
collection.  Item  life  is  measured  in  terms  of  cycles  of  operation.  Protection 


againdt  lots  for  which  the  average  item  life  is  less  thar.  5>000  cycles  is 
required.  From  experience  with  this  product  it  has  been  determined  that  the 
Weibull  distribution  applies  and  that  a  value  for  7  of  2,000  cycles  and  a  value 
for  3  or  2  can  be  expected.  The  problem  is  to  determine  a  teat  time,  t,  in 
cycles  that  will  enable  the  plan  to  meet  the  above  requirement  for  consumer's 
protection.  A  related  problem  is  to  find  whether  the  plan  so  determined  will  give 
adeouate  producer  protection.  It  has  been  determined  that  a  mean  item  life  of 
16,000  can  reasonably  be  expected  from  a  competent  supplier . 

The  first  step  toward  a  solution  is  to  convert  the  required  lot  mean  life, 
u  to  a  transformed  value,  This  new  value,  is  u  -  7  or  000  -  2,000 

which  is  3,000  cycles.  Next,  from  Table  5f  which  gives  conversion  ratios  for  use 
when  3=2,  one  finds  that  for  Sample  Size  Letter  H  and  an  AQL  in  p'  (^)  of  l-5> 
the  ratio  at  the  LTPD  Quality  Level  is  38  end  at  the  AQL  it  is  13-87-  Since  the 
plan  is  to  be  determined  in  terms  of  consumer's  needs,  the  next  step  is  to  use 

the  LTPD  ratio  to  determine  t  .  Thus  (t  /u  )  x  100  ■  38  or  (t  /3>°00)  x  100  »  38. 

000  y 

From  this  it  is  determined  that  t  must  equal  1140.  3y  adding  the  value  for  7 
(which  is  2,000)  to  this  latter  figure,  the  required  test  time  in  absolute  terms, 
t,  of  3,140  cycles  is  obtained. 

The  related  question  of  the  reasonableness  of  this  plan  for  the  producer 
may  be  answered  by  substitution  of  the  test  time  Just  determined  in  the  conversion 
ratio  for  the  AQL.  The  relationship  is  that  (t0/|ie)  x  100  *  13-87  or  that 
(114o/uq)  x  100  -  13.87.  From  this  a  value  for  uQ  of  8,240  cycles  is  obtained. 
This  is  then  converted  to  original  terms  by  the  addition  of  the  value  for  7  of 
2.000  cycles.  This  gives  an  Acceptable  Quality  I^vel  of  10,240  cycles.  This  is 
well  below  the  level  considered  reasonable  so  no  hardship  will  be  impojed  on  the 
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TABLE  1 

Table  of  Values  for  Per  cent  Truncation,  (t/n)  x  100 


p’ 

(in  %) 

Shape  Parameter  -  8 

V3 

1/2 

1 

1  2/3 

2 

2  1/2 

3  1/3 

U 

.010 

.010 

•U5 

1.13 

2.93 

7.03 

11.03 

17.26 

.012 

.012 

.L9 

1.2U 

3-OU 

7.U2 

11.55 

17.91 

.015 

.015 

.57 

1.36 

3.32 

7.9U 

12.21 

18.72  1 

.020 

.020 

.67 

1.59 

3.73 

8.66 

13.12 

19.83 

.025 

.025 

.77 

1.78 

U.08 

9.26 

13.87 

20. 7U 

.030 

.030 

.96 

1.95 

U.UO 

9.77 

1U.52 

21.50 

.01*0 

.ouo 

1.02 

2.26 

U.93 

10.65 

15.60 

22.77 

.050 

.050 

1.19 

2.53 

5.39 

11. UO 

16.  U9 

23.82  | 

.065 

.065 

1.37 

2.88 

5.98 

12.32 

17.62 

25.10  ! 

.080 

.080 

1.56 

3.19 

6.50 

13.13 

18.56 

26.16  | 

.100 

.10 

1.78 

3.57 

7.11 

1U.03 

19.62 

27.36  ! 

.12 

.12 

1.98 

3.92 

7,65 

1U.32 

20.53 

28.37 

.15 

.15 

2.26 

U.37 

8.36 

15. 8U 

21.71 

29.67 

.20 

.20 

2;69 

5.07 

9.39 

17-27 

23.33 

31.  U3 

.25 

.25 

3.08 

5.6U 

10.27 

18. U7 

2U.68 

32.87 

.30 

.30 

3.UU 

6.18 

11.05 

19.51 

25.83 

3U.09 

.Uo 

.uo 

U.07 

7.1U 

12.39 

21.27 

27.76 

36.12 

.50 

.001 

.50 

U.67 

7.99 

13.55 

22.75 

29.36 

37.76 

.6 5 

.002 

.65 

5.U6 

9.12 

15.06 

2U.62 

31.35 

39.81 

.80 

.003 

.80 

6.19 

10.11 

16.36 

26.21 

33.03 

U1.50 

1.00 

.005 

1.01 

7.08 

11.31 

17.90 

28.03 

3U.93 

U3.U0 

1.2 

.007 

1.21 

7.90 

12. UO 

19.26 

29.62 

36.57 

U5.02 

1.5 

.on 

1.51 

9.07 

13.87 

21.08 

31.68 

38.68 

U7.09 

2.0 

.020 

2.02 

10.77 

16.03 

23.67 

3U.56 

U1.59 

U9.90 

2-5 

.032 

2.53 

12.33 

17.95 

25.90 

36.98 

UU.01 

52.21 

3.0 

.0U7 

3.05 

13.78 

19.69 

27.89 

39.09 

U6.09 

5U.17 

L.O 

.001 

.083 

u.08 

16.  U2 

22.79 

31.35 

U2.69 

U9.59 

57.U5 

5.0 

.002 

.13 

5.13 

18. 8U 

25.58 

3U.35 

U5.71 

52.50 

60.13 

6.5 

.005 

.23 

6.72 

22.15 

29.25 

38.28 

U9.57 

56.18 

63.U6 

8.0 

.010 

.35 

8.3U 

25.20 

32.59 

Ul.72 

52.88 

59.29 

66.26 

10.0 

.020 

.56 

10. 5U 

29.01 

36.63 

U5.82 

56.73 

62.85 

69.UU 

12 

.03U 

.82 

12.78 

32.58 

U0.3U 

U9.50 

60. 11 

65.96 

72.18 

15 

.070 

1.32 

16.25 

37.63 

U5.U6 

5U.U9 

6U.6O 

70.05 

75.73 

20 

.18 

2.U9 

22.31 

U5.51 

53.30 

61.85 

71. OU 

75.83 

80.68 

25 

.Uo 

U.1U 

28.77 

52.99 

60.53 

68. U7 

76.67 

80.80 

8it.89 

30 

.76 

6.36 

35.37 

60.29 

67.39 

7U.62 

81.79 

85.26 

88.62 

uo 

2.22 

13.  OU 

51.08 

7U.79 

80. 6U 

86.15 

91.09 

93.27 

95.22 

50 

5.55 

2U.02 

69.31 

89.92 

93.95 

97.33 

99.82 

100.67 

101.21 

65 

19.28 

55.10 

10U.98 

115*23 

115. 61 

11U.92 

113.06 

111.68 

109.98 

80 

69.U8 

129.52 

160. 9U 

1U8.91 

1U3.1U 

136.3U 

128.53 

12U.27 

119.79 
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TABLE  2 


Table  of  Probability  Values  at  Truncation  Point,  p'  (%) 


Shape  Parameter  -  0 


(t/n)  X  100 

V3 

1/2 

1 

1  2/3 

2 

2  1/2 

.010 

>- 

8.09 

1.1*0 

.010 

.012 

8.57 

1.51* 

.012 

.015 

9.20 

1.72 

.015 

.020 

!  io.o8 

1.98 

.020 

.025 

1  10.82 

2.21 

.025 

.030 

11.1*5 

2.1*2 

.030 

.01*0 

12.53 

2.79 

.01*0 

.050 

13.1*3 

3.11 

.050 

.065 

11*.56 

3.51* 

.065 

.080 

15.52 

3.92 

.000 

.100 

16.61 

U.37 

.10 

.12 

17.56 

1*.78 

.12 

.001 

.15 

18.78 

5.33 

.15 

.002 

.20 

20.1*6 

6.13 

.20 

.003 

.25 

21.86 

6.83 

.25 

.001* 

.30 

23.06 

7.1*5 

.30 

.005 

.1*0 

25.06 

8.56 

•1*0 

.009 

.001 

.50 

26.71 

9.52 

.50 

.012 

.002 

.65 

28.76 

10.78 

.65 

.019 

.003 

.80 

30.1*7 

11.88 

.80 

.027 

.005 

1.00 

32. 1*0 

13.19 

1.00 

.038 

.008 

1.2 

31*. 03 

11*. 35 

1.19 

.052 

.on 

.001 

1.5 

36.12 

15.90 

1.1*9 

.076 

.018 

.002 

2.0 

38.91* 

18.13 

1.98 

.12 

.031 

.001* 

2.5 

U1.22 

20.01* 

2.1*7 

.18 

.01*9 

.007 

3.0 

1*3.11* 

21.73 

2.96 

.21* 

.071 

.012 

L.o 

1*6.28 

21*. 61* 

3.92 

.39 

.13 

.021* 

5.0 

1*8.80 

27.11 

1*.88 

.56 

.20 

.ota 

6.5 

50.90 

30.27 

6.29 

.89 

.33 

.080 

8.0 

51*.  30 

32.97 

7.69 

1.22 

.50 

.13 

10.0 

56.98 

36.06 

9.52 

1.77 

.78 

.23 

12 

59.19 

38.73 

11.31 

2.39 

1.12 

.37 

15 

61.92 

1*2.17 

13.93 

3.1*5 

1.75 

.61* 

20 

65.1*5 

1*6.87 

18.13 

5.51 

3.09 

1.32 

25 

68.17 

50.69 

22.12 

7.89 

1*.79 

2.29 

30 

70.37 

53.91 

25.92 

10.56 

6.82 

3.59 

1*0 

73.79 

59.12 

32.97 

lc.l*7 

11.81 

7.23 

50 

76.36 

63.21 

39.35 

22.98 

17.83 

11.28 

65 

79.28 

68.02 

1*7.80 

33.26 

28.21* 

22.3  2 

80 

81.1*9 

71.77 

55.07 

1*3.51* 

39.51 

31*. 59 

100 

|  83.75 

75.69 

63.21 

56.35 

51*. la 

52.36 

15.37  11.35 

28.35  21*.  15 

50.1a  1*9.08 


2. 
7.29 
19.25 
1*7.93 
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TABLE  3a 

Tabla  of  Stapling  Plana  for  6  •  1/3 


(t/**)  x  100  Ratio  for  which  P(A)  •  *10  (or  Leas) 


25  10 


3  1*  h  6 

(.05)  (.02)  (.01) 

*  5  6  8 

(.15)  (.16)  (.08)  (.03) 

6  7  8  10 

(.53)  (.27)  (.17)  (.07) 

7  8  10  12 

(1.2)  (.66)  (.33)  (.Hi) 

9  10  11  U* 

(-1.3)  (.80)  (.51*)  (.20) 


0.5 

0.25 

0.1 

0.05 

10 

(2.2) 
11 

(2.6) 

13 
(2.9) 

9  11* 

(3.8) 

15 

(1*.6) 


11  16 
(5.0) 

12  18 
(5.5) 

13  19 

(6.2) 

11*  20 
(6.7) 

15  22 

(7.0) 


11  13 

(1.2)  (.65) 

13  15 

(1.1*)  (.76) 

11*  16 
(1.7)  (1.1) 

16  18 

(2.1)  (1.2) 

17  20 
(2.5)  (1.2) 


19 

(2.8) 

20 

(3.1)  . 

21  21* 
(3.6)  (1.9) 

23  26 
(3.9)  (2.1) 

21*  28 

(U.2)  (2.2) 


26 

(.68) 

28 

(.72) 

30 

(.76) 

32 

(.85) 

31* 

(.95) 


(.23) 

27 

(.21*) 

29 

(.28) 


31 

(.31*) 

31* 

(.31*) 

36 

(.1*0) 

38 

(.1*5) 

1*0 

(.1*5) 


31  37  1*6 

(.03)  (.02)  (.01) 

31*  1*2  52  oy 

(.01*)  (.02)  (.01) 

38  1*7  58  76 

(.05)  (.02)  (.01) 

1*1  51  63  83 

(.06)  (.03)  (.01)  (.01) 

1*5  55  71  93 

(.06)  (.03)  (.01)  (.01) 


60 

(.03) 


(t/}i)  x  100  ratios  in  paxanthasaa  art  for  P(A)  ■  .95  (or  aora) 
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TABLE  3b 


Table  of  Sampling  Plana  for  8  ■  X/2 


1*77 

753 

533 

81*1 

589 

929 

61*3 

1,015 

698 

1,101 

752 

1,186 

805 

1,271 

858 

1,355 

911 

(.01) 

1,1*38 

961* 

(.01) 

1,5a 

(t/p)  x  100  ratios  in  parentheses  are  for  P(A)  •  .95  (or  more) 
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TABLE  3c 


r 

Tabla  of  Sampling  Plan®  for  £  •  1 


(t/ti)  z  100  ratio®  in  p®r«nth®®es  ar®  for  P(A)  •  .95  (or  aor®) 

Th®  figur®  following  th®  dash  in  sample  sis®  number®  show®  th®  nuitbar  of  s®ro®  to 
add;  for  ®xaapl®j  15L-3  •  15U,000. 
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TABLE  3d 

Table  of  Sailing  Plana  for  0  ■  1  2/3 


(t /n)  x  100  Ratio  for  which  P(A)  ■  .10  (or  laaa) 


a 

15 

■ 

8 

_ * 

5 

■3 

2 

s 

■ 

B 

5  o.25 


11 

26 

(59) 

12 

28 

(60) 

13 

30 

(61) 

u* 

33 

(61) 

15 

35 

(62) 

1*37-2  133-3  105-3 
(.53)  (.30)  (05) 

1*66-2  11*2-3  1*15-3 
(.59)  (.30)  (.15) 

1*99-2  152-3  1*75- 
(.61)  (.31)  (.15 

530-2  161-3  505-3 
(.62)  (.31)  (.16 

116-2  260-2  560-2  170-3  535- 
(.16)  (.95)  (.63)  (.32)  (.16 


(t/t*)  x  100  ratioa  in  paranthaaas  ara  for  P(A)  ■  .95  (or  more) 

Tha  figure  following  the  daah  in  staple  aiaa  numbers  shows  the  number  of  Mros 
to  add}  for  example,  216-2  •  21,800. 


2 
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TABLE  3e 

Table  of  Sampling  Plana  for  6  ■  2 


n 

0 

— 

(t/**) 

x  100  Ratio  for  which  P(A)  -  * 

>10  (or  less) 

100 

50 

25 

15 

12 

10 

8 

5 

1* 

2.5 

1.5 

1 

0.5 

0 

3 

12 

1*7 

131 

206 

296 

1*61 

1,152 

1,772 

1*,700 

128-2 

208-2 

11^-3 

kis) 

7.U) 

(3.7) 

(2.2) 

(1.8) 

(1.5) 

|(1.2) 

(.76) 

(.62) 

(.38) 

(.23) 

(.16) 

(.08) 

l 

6 

21 

80 

223 

31*6 

1*99 

778 

1,91*5 

2,993 

7,939 

216-2 

1*86-2 

19l*-3 

<29) 

(15) 

(7.5) 

(1*.5) 

(3.6) 

(3.0) 

(2.1*) 

(1.5) 

(1.2) 

(.76) 

(.1*6) 

(.30) 

(.15) 

2 

8 

29 

no 

305 

1*76 

683 

1,061* 

2,661 

1*,091* 

109-2 

296-2 

665-2 

266-3 

k38) 

(19) 

(9.8) 

(5.8) 

(i*.7) 

(3.9) 

(3.1) 

(2.0) 

(1.6) 

(.98) 

(.60) 

(.1*0) 

(.20) 

3 

11 

36 

138 

382 

597 

857 

1,336 

3,31*1 

5,H*o 

136-2 

371-2 | 835-2 

331*-3 

(1*3) 

(22) 

(U) 

(6.7) 

(5.1*) 

(1*.  5) 

(3.6) 

(2.3) 

(1.8) 

(1.1) 

(.70) 

(.1*7) 

(.21*) 

1* 

13 

1*3 

167 

1*57 

711* 

1,025 

1,599 

3,977 

6,150 

163-2 

1*1*1*— 2  999-2 

1*00-3 

(U8) 

(21*) 

(12) 

(7.1*) 

(5*9) 

(5.0) 

(H.O) 

(2.5) 

(2.0) 

(1.2) 

(.76) 

(.51) 

(.26) 

5 

15 

50 

191* 

530 

829 

1,190 

1,855 

1*,638 

7,135 

189-2 

515-2  116-3 

l*61*-3 

(52)|  (26) 

(13) 

(7.9) 

(6.3) 

(5.3) 

(1*.2) 

(2.7) 

(2.2) 

(1.3) 

(.81)  (.51*)  '(.28) 

6 

1  17 

57 

220 

602 

91*1 

1,351 

2,106 

5,211* 

8,102 

215-2 

585-2! 132-3 

527-3 

(55) 

(28) 

(11*) 

(8.3) 

(6.6) 

(5.6) 

(1*.5) 

(2.8) 

(2.3) 

(1.1*) 

(.85) 

(.57) 

(.29) 

7 

19 

61* 

21*5 

673 

1,051 

1,510 

2,351* 

5,806 

9,055 

21*0-2 

65L-2  11*7-3 

589-3 

(58) 

(29) 

(U*) 

(8.7) 

(6.9) 

(5.8) 

(1*.6) 

(2.9) 

(2.1*) 

(1.5) 

(.88) 

(.60) 

(.30) 

8 

21 

71 

272 

71*3 

1,161 

1,667 

2,600 

6,1*98 

9,997 

265-2 

722-2 

162-3 

650-3 

(60) 

(30) 

(15) 

(9.0) 

(7.2) 

(6.0) 

(1*.8) 

(3.0) 

(2.1*) 

(1.5) 

(.91) 

(.62) 

(.31) 

9 

23 

77 

297 

812 

1,271* 

1,822 

2,81*1 

7,103 

109-2 

290-2 

789-2 

178-3 

710-3 

(62) 

(30) 

(15) 

(9.2) 

(7.1*) 

(6.1) 

(1*.9) 

(3.1) 

(2.5) 

(1.5) 

(.91*) 

(.63) 

(.32) 

10 

26 

87 

322 

881 

1,376 

1,976 

3,081 

7,701* 

119-2 

3H*-2 

856-2 

193-3 

770-3 

(62) 

(31) 

(15) 

(9.1*) 

(7.5) 

(6.3) 

(5.0) 

(3.2) 

(2.6) 

(1.6) 

(.96) 

(.65) 

(.33) 

11 

28 

91* 

31*7 

91*9 

1,1*82 

2,128 

3,320 

8,299 

126-2 

339-2 

922-2 

20703 

830-3 

(63) 

(31) 

(16) 

(9.6) 

(7.7) 

(6.1*) 

(5.1) 

(3.3) 

(2.6) 

(1.6) 

(9.8) 

(.66) 

(.31*) 

12 

30 

100 

372 

b.,017 

1,588 

2,280 

3,556 

8,891 

137-2 

363-2,988-2 

222-3 

009-3 

(61*) 

(32) 

(16) 

(9.8) 

(7.8) 

(6.5) 

(5.2) 

(3.3) 

(2.7) 

(1.6) 

(1.0) 

(.67) 

(.31*) 

13 

32 

107 

396 

1,081* 

1,693 

2,1*31 

3,792 

9,1*79 

11*6-2 

387-2 

105-3 

237-3  '91*8-3 

(65) 

(32) 

(16) 

t.  do) 

(8.0) 

(6.6) 

(5.3) 

(3.L) 

(2.7) 

(1.7) 

(1.0) 

(.68)!  (.35) 

lit 

31* 

113 

1*21 

1,151 

1,798 

2,581 

1*,026 

101-2 

155-2 

1*11-2 

112-3 

252-3 

101-1* 

(66) 

(33) 

(17) 

(10) 

(8.1) 

(6.7) 

(5.1*) 

(3.1*) 

(2.8) 

(1.7) 

(1.0) 

(.69) 

(.35) 

15 

36 

120 

1*1*5 

1,217 

1.902 

2,730 

1*,258 

106-2 

l61*-2 

1*35-2 

118-3 

266-3 

106-1* 

(67) 

(33) 

(17) 

(10) 

(8.2) 

(6.8) 

(5.5) 

(3.5) 

(2.8) 

(1.7) 

(1.0) 

(.70) 

(.36) 

(t/**)  x  100  ratios  in  partnthases  are  for  P(A)  •  .95  (or  more) 


The  figure  following  the  daeh  In  sample  size  numbers  shows  ti>e  number  of  zeros 
to  add)  for  example,  118-3  •  118,000* 
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table  3f 

Table  of  Sampling  Plans  for  £  •  2  1/2 


(t/p)  x  100  Ratio  for  which  P(A)  ■  *10  (or  lass) 


320-2  115-3 
(.55)  (.33) 

51*0-2  19l*-3 
(.96)  (.58) 

739-2  266-3 
(1.2)  (.71) 


51*1-2 

(2.7) 

592-2 

(2.7) 

61*2-2 

(2.8) 


(1.3) 

111-3 

(1.1*) 

129-3 

(1.5) 


527-3 
(.91*) 

589-3 
(.97) 

180-3  650-3 
(1.6)  (1.0) 

197-3  710-3 
(1.7)  (1.0) 

2H*-3  770-3 
(1.7)  (1.0) 


692-2  231-3 

(2.8)  (1.7) 

71*1-2  21*7-3 

(2.9)  (1.8) 

790-2  263-3  91*8-3 
(2.9)  (1.8)  (1.1) 

839-2  280-3  101-1* 
(2.9)  (1.8)  (1.1) 

887-2  296-3  106-1* 
(3.0)  (1.8)  (1.1) 


(t/n)  x  100  ratios  in  parentheses  ars  for  P(A)  ■  .95  (or  mora) 

Ths  figura  following  the  dash  in  sampla  sisa  numbs rs  shows  tha  number  of 
seros  to  add)  for  example,  296-3  ■  296,000. 
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TABLo  3g 

Table  of  Sampling  Plana  for  0  ■  3  1/3 


■‘"/u)  x  100  ratio*  in  parent  base*  are  for  P(A)  •  • 9 $  (or  more) 

The  figure  following  the  daeh  in  sample  site  number*  shove  the  nvaber  of 
taro*  to  add)  for  example,  319-3  •  319,000. 
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TABLE  3h 

Table  of  Sampling  Plans  for  P  ■  U 


n 

(t/*0  X 

LOO  Ratio  for  which  P(A) 

-  .10 

(or  less) 

100 

T 

80 

65 

50 

UO 

P 

— 

25 

20 

l 

15 

12 

i 

10 

8 

1  6.5 

j  _ . 

0 

U 

9 

20 

■p” 

55 

13U 

U19 

886 

2,09U 

,  6,77U 

16U-2 

329-2 

768-2 

177-3 

(37) 

(30) 

1  (25) 

1  (19) 

(15) 

i  (12) 

(9.U) 

(7.6) 

(5.8) 

(U.7) 

(U.o) 

(3.2) 

(2.6) 

1 

7 

15 

1 

33 

93 

228 

7oe 

1,U97 

3,537 

HU-2 

278-2 

556-2 

130-3 

299-3 

(53) 

(UU) 

(36) 

;  (27) 

(22) 

(16) 

| 

(1U) 

(11) 

(8.0) 

(6.5) 

(5.6) 

(U.5) 

(3.7) 

2 

9 

21 

1*6 

j  128 

312 

968 

2,OU7 

U,839 

157-2 

380-2 

760-2 

177-3 

U09-3 

(62) 

(50) 

(UO) 

|  (31) 

(25) 

(19) 

(16) 

'  (12) 

(9.2) 

(7.U) 

(6.2) 

(5.1) 

(U-2) 

3 

12 

26 

57 

162 

391 

1 1,215 

2,570 

6,07U 

197-2 

U77-2 

95U-2 

223-3 

51U-3 

(66) 

<5U) 

m 

(33) 

;  (27) 

(20^ 

(17) 

(13) 

(10) 

(7.9) 

(6.6) 

(5.5) 

(U.5) 

U 

15 

31 

69 

;  191* 

U66 

l,U5U 

3,075 

7,268 

235-2 

571-2 

llU-3 

266-3 

1 615-3 

(68) 

(56) 

<U6) 

(35) 

(28) 

(21) 

(18) 

(1U) 

(10) 

(8.3) 

(7.0) 

(5.8) 

(U.7) 

5 

17 

37 

00 

225 

5U3 

1,687 

3,568 

8,1*32 

273-2 

663-2 

133-3 

309-3 

713-3 

(71) 

(58) 

(1*8) 

(36) 

(29) 

(22) 

(18) 

(15) 

(11) 

(8.6) 

(7.2) 

(6.0) 

(U.9) 

6 

19 

U 2 

91 

255 

616 

1.915 

U,051 

9,575 

310-2 

t 

752-2 

150-3 

351-3 

810-3 

(7U) 

(59) 

(1*9) 

(37) 

(30) 

(22) 

(19) 

(15) 

(11) 

(8.8) 

(7.U) 

(6.1) 

(5.0) 

7 

22 

U7 

102 

286 

689 

2,1U0 

U,520 

107-2, 

3U6-2 

81*1-2 

168-3 

392-3 

905-3 

(75) 

(60) 

(50) 

(38) 

(30) 

(23) 

(19) 

(15) 

(11) 

(9.0) 

(7.6) 

(6.2) 

(5.1) 

0 

2U 

52 

112 

315 

760 

2,363 

U,998 

118-2 

382-2 

928-2 

186-3 

U33-3 

100-U 

(76) 

(61) 

(50) 

(39) 

(31) 

(23) 

(19) 

(15) 

(12) 

(9.1) 

(7.7) 

(6.3) 

(5.2) 

9 

27 

57 

123 

31*1* 

831 

2,583 

5,U6U 

129-2 

1*18-2 

101-3 

203-3 

U7U-3 

109-u 

(77) 

(62) 

(51) 

(39) 

(31) 

(23) 

(20) 

(16) 

(12) 

(9.2) 

(7.8) 

(6.3) 

(5.3) 

10 

29 

61* 

136 

373 

901 

2,802 

5,926 

11*0-2 

U53-2 

110-3 

220-3 

51U-3 

119-U 

(78) 

(62) 

(51) 

(1*0) 

(32) 

(2U) 

(20) 

(16) 

(12) 

(9.U) | 

(7.8) 

(6.U) 

(5.3) 

11 

31 

69 

11*7 

1*02 

971 

3,018 

6,38U 

151-2 

U88-2 

'T 

119-3 

237-3 

553-3 

128-1* 

(79) 

(63) 

(52) 

(Uo) 

(32) 

(2U) 

(20) 

(16) 

(12) 

(9.U)1 

(7.9) 

(6.U) 

(5.U) 

12 

33 

71* 

157 

1*31 

i,ouo 

3,233 

6,8U0 

162-2 

523-2 

127o! 

25U-3 

593-3 

137-U 

(80) 

(63) 

(52) 

(UO) 

(32) 

(2U) 

(20) 

(16) 

(12) 

(9.5) 

(8.0) 

(6.5) 

(5.U) 

13 

36 

79 

167 

1*60 

1,109 

3,UU7 

7,292 

172-2 

558-2 

135-3 1 

271-3 

632-3 

1U6-U 

(81) 

(61*) 

(53) 

(33) 

(2U) 

(20) 

(12) 

(9.6)! 

(8.1) 

(6.5) 

(5.5) 

iu 

38 

81* 

178 

U88 

1,177 

3,660 

7,7U2 

193-2 

592-2 

1UU-3, 

286-3 

671-3 

155-U 

(82) 

(61*) 

(53) 

(U) 

?33) 

(25) 

(21) 

(17) 

(12) 

(9*7) | 

(8.2) 

(6.6) 

(5.5) 

15 

L_ 

Uo 

89 

188 

516 

1,21*6 

3,872 

8,190 

19U-2 

626-2 

152-3 

30U-3 

710-3 

16U-U 

(82) 

(61*) 

(53) 

(Ul) 

(33) 

(25) 

(21) 

(17) 

(12) 

(9.8)1 

(8.2) 

(6.6) 

(5.6) 

(t/V)  x  100  ratios  In  parentheses  are  for  P(A)  »  .95  (or  more) 

The  figure  following  the  dash  in  sample  site  nunbera  shows  the  mmber  of 
aeros  to  add;  for  example,  30U-3  •  30U,000. 
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TARlx  3« 

Table  of  Sampling  Plans  for  (3  ■  3  1/3 


n 

(t/t*)  x  100  Ratio  for  which  P(a) 

•  .10  (or  less) 

100 

1  65 

50 

1*0 

3° 

25 

1  _  .  _ 

20 

15 

12 

10 

8 

6.5 

5 

0 

1* 

11* 

31* 

70 

! 

|  183 

331* 

698 

‘1,772 

3,839 

6,979 

11*9-2 

288-2 

698-2 

(30) 

!  (21) 

(16) 

(13) 

(9.6) 

(8.0) 

(6.1*) 

(1*.8) 

(3.8) 

|  (3*2) 

(2.5) 

(2.1) 

(1.6) 

1 

7 

21* 

57 

119 

309 

561* 

1,179 

2,993 

6,1*81* 

!  118-2 

251-2 

1*86-2 

118-3 

(1*6) 

I  (32) 

(21*) 

(19) 

I  (H*) 

(12) 

(9.8) 

1  (7.1*) 

(5.8) 

kii.s) 

(3.9) 

(3.2) 

(2.1*) 

2 

9 

33 

78 

161* 

1  1*23 

772 

1,613 

l*,09l* 

8,870 

1 161-2 

31*3-2 

665-2 

161-3 

(56) 

(37) 

(28) 

1  (23)  I  (17) 

(H*) 

(11) 

(8.7) 

(6.8) 

j (5*7) 

(1**5) 

(3.7) 

(2.8) 

3 

12 

1*2 

98 

206 

531 

969 

2,025 

5,11*0 

111-2 

202-2 

1*31-2 

1835-2 

202-3 

(60) 

(1*0) 

(31) 

(25) 

(19) 

(15) 

(12) 

(9.1*) 

(7.5) 

(6.2) 

(1**9) 

,(1*.0) 

(3.1) 

1* 

li* 

51 

117 

21*6 

635 l 1,159 

2,1*23 

6,150 

133-2 

21*2-2 

516-2 

999-2 

21*2-3 

i 

is 

1  _ 

(65) 

(1*2) 

(33) 

(26) 

(20) 

(16) 

(13) 

(10) 

(7.9) 

(6.6) 

(5.2) 

(1*.3) 

(3.3) 

16 

(68) 

59 

(U*) 

136 

(31*) 

266 

1(27) 

737 

(21) 

1,31*5 

(17) 

2,811 1 7.135 
(H*)  1  (10) 

155-2 

(8.2) 

281-2 

(6.8) 

598-2 

(5.1*) 

116-3 

(1*.U) 

281-3 

(3.U) 

» 

6 

19 

67 

157 

325 

836 

1,527 

3,192 

8,102 

176-2 

319-2 

679-2 

132-3 

319-3 

(69) 

(1*6) 

(35) 

(28) 

(21) 

(17) 

(U*) 

(11) 

(8.1*) 

(7.0) 

(5.6) 

(1*.6) 

(3.5) 

7 

21 

75 

176 

363 

935 

1,698 

3,567 

9,055 

196-2 

357-2 

759-2 

11*7-3 

357-3 

(71) 

(1*7) 

(36) 

(29) 

(22) 

(18)  (11*) 

(11) 

(8.7) 

(7.2) 

(5.8) 

(U.7) 

(3.6) 

e 

23 

83 

191* 

1*00 

1,032 

1,881* 

3,93b 

9,997 

217-2 

39U-2 

b3b-2 

lfec-3 

391*-3 

(73) 

(1*7) 

(36) 

(29) 

(22) 

(18) 

(15) 

(U) 

(8.8) 

(7.1*) 

(5.9) 

(1*.8) 

(3.7) 

9 

26 

90 

212 

1*38 

1,128 

2,059 

1*,305 

109-2 

237-2 

1*30-2 

917-2 

178-3 

1*30-3 

(7U) 

(1*8) 

(37) 

(30) 

(22) 

(19) 

(15) 

(11) 

(9.0) 

(7.5) 

(6.0) 

(1**9) 

(3.7) 

10 

26 

101 

230 

1*75 

1,228 

2,233 

1*,669 

119-2 

257-2 

1*67  -2 

99l*-2 

193-3 

1*67-3 

(75) 

(1*8) 

(38) 

(30) 

(23) 

(19) 

(15)1 

(12) 

(9.2) 

(7.6) 

(6.0) 

(1*.9) 

(3.8) j 

li 

30 

108 

21*7 

511 

1,318 

2,1*06 

5,030 

128-2 

277-2 

503-2 

107-3 

207-3 

503-3 

(76) 

(1*9) 

(38) 

(31) 

(23) 

(19) 

(15) 

(12) 

(9.3) 

(7.7) 

(6.1) 

(5.0) 

(3.8) 

12 

32 

116 

265 

51*8 

1,1*12 

2,578 

5,389 

137-2 

296-2 

539-2 

115-3 

222-3 

539-3 

(77) 

(1*9) 

(39) 

(31) 

(23) 

(20) 

(15) 

(12) 

(9.1*) 

(7.8) 

(6.2) 

(5.0) 

(3.9) 

13 

35 

121* 

283 

581*  1,505 

2,71*8 

5,71*5 

11*6-2 

316-2 

57U-2 

122-3 

237-3 

571*— 3 

(77) 

(50) 

(39) 

(31)| 

(21*)  ( 

(20) 

(16) 

(12) 

(9.U) 

(7.9) 

(6.2) 

(5.1) 

(3.9) 

11* 

37 

131 

300 

620 

1*596 

2,918 

6,100 

155-2 

335-2 

610-2 

130-3 

252-3 

610-3 

(78) 

(50) 

(39) 

(32) 

(21*) 

(20)i 

(16) 

(12) 

(9.5) 

(8.0) 

(6.3) 

(5.2) 

(3.9) 

15 

39 

139 

317 

656 

1,690 

3,086;  6,1*53 

l61*-2 

355-2 

61*5-2 

137-3 

266-3 

61*5-3 

(79) 

(51) 

(UO) 

(32) 

(21*) 

(20)| 

_ L 

(16) 

(12) 

(9.6) 

(8.0) 

(6.1*) 

(5.2) 
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Tha  figure  following  the  dash  in  sample  site  numbers  shows  the  number  of 
taros  to  add;  for  example,  319-3  •  319,000. 
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The  figure  following  the  dash  in  sample  site  numbers  shows  ths  numbs r  of 
seros  to  add;  for  example ,  30L-3  ■  30L,000. 
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TABLE  31 

Table  of  Samplir^  Plane  for  P  •  5 


(t/V)  x  100  ratios  in  parsnthesss  ars  for  P(A)  •  ,95  (or  more) 

The  figure  following  the  dash  in  sample  size  numbers  shows  the  number  of 
zeros  to  add}  for  example,  203-2  ■  20,300. 
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Table  of  (t/ti)xlOO  Ratios  for  0  -  l/2  for  the  105B  Plans 
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Table  of  (t/ii)xl0O  Ratios  for  p  =  2  for  the  105B  Plans 
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APPENDIX  A 


the  Welbull  Distribution  as  a  Failure  Distribution 

This  appendix  describes  the  Weibuil  distribution  which  is  used  as  a 
mathematical  model  for  the  procedures  described  in  this  report.  The  Wv'ibull 
distribution  is  a  member  of  the  general  failure  distribution  which  can  be 
derived  as  follows: 

Let  N  »  the  initial  number  of  items  under  life  testing. 

S(x)  =  the  number  of  items  still  surviving  at  time  x. 

R(x)  =  S(x)/N  =  the  proportion  of  items  still  surviving  at  time  x. 
F(x)  *  1  -  R(x)  -  the  proportion  of  failures  prior  to  time  x. 

Z(x)  =  the  Instantaneous  rate  of  failure  at  time  x. 

7  -  An  assumed  time,  prior  to  which  no  failure  will  occur. 

The  introduction  of  7,  some  finite  time,  means  when  x  t  7 ;  S(x)  =  N,  R(x)  ■  1, 
F(x)  =  Z(x)  *  0  which  is  a  boundary  condition.  A  negative  value  for  7  means 
that  the  items  could  fall  prior  to  life-testing,  i.e.,  the  items  could  fail  in 
storage,  a  common  phenomenon  for  batteries,  electrolytic  capacitors,  electron 
tubes  and  the  like. 

The  following  simple  differential  equation  equates  the  decrement  of  items 
surviving  in  some  small  increment  of  time  dx  with  the  product  of  S(x)  and  Z(x). 
The  negative  sign  indicates  the  decrement. 

dS(x)/dx  -  -Z(x)  •  S(x)  (A  l) 

Divide  both  sides  by  n  and  transpose, 

dR(x)/R(x)  -  -Z( x)  dx  (A  2) 

Now  integrate  both  sides  over  the  interval  (  7,  x  ), 

X  X 

In  R(y)  I  -  -  /  Z(y)  dy  +  C  (A3) 

7  7 

The  integration  constant  C  is  zero,  because  of  the  boundary  condition  of 
R(x)  =1  and  Z(x)  =  0,  for  x  S  7 
Hence , 

x 

In  R(x)  -  -  /  Z(y)  dy  (A  4) 

7 

-Vfr 


which  can  be  rewritten  as 

x 

F(x)  «  1  -R(x)  -  1  -  exp  (  -  /  Z(y)dy  ]  »  1  -  exp  [  -M(x)  1  f or  x  i  7, 

7 

*  0  ,  for  x  <  7  .  (A  5) 

Ibis  is  the  form  that  a  failure  distribution  must  take.  Since  Z(x)  is  non* 

negative  for  all  x,  in  order  lor  F(x)  to  have  meaning,  M(x)  oust  be  non- 

decreasing  in  x.  Although  there  are  still  infinitely  many  choices  of  M(x), 

this  form  of  distribution  rules  out  many  other  distributions,  for  example  the 

normal  distribution,  as  a  failure  distribution. 

If  the  Instantaneous  failure  rate,  Z(x)  is  a  constant  over  time  say  equal 
x 

to  X,  then  M(x)  *  |  X'dy  «  X(x-7),  for  x  >  7,  and  M(x)  ■  0  for  x  S  7. 

7 

The  resulting  distribution  is  the  well-known  2 -parameter  exponential  dlstributiot 
namely  for  X  ■  1/0, 

F(x)  ■  1  -  exp  [  -(x-7)/0  ]  ,  for  x  *  7,  $  >  0, 

■  0,  otherwise.  (A  6) 

However,  if  the  instantaneous  failure  rate  changes  with  time  in  such  a  way  that 
.y\P 

M(x)  «  — for  x  >  7  and  M(x)  -  0  for  x  S  7,  the  resulting  distribution 
will  be  the  following: 

F(x)  *  1  -  exp  (  -(x-7)^/a  ]  ,  for  x  ft  7,  a,  p  >  0, 

■  0,  otherwise.  (A  7) 

This  form  of  the  failure  distribution  is  known  as  the  Welbull  distri¬ 
bution.  The  three  parameters  of  the  Welbull  distribution  are, 
a  -  the  "scale"  parameter 
0  •  the  shape  parameter 

7  -  the  location  parameter 

Often  the  Welbull  distribution  is  written  as, 

P  * 

F(x)  -  1  -  exp  [  -(-£lL)  ]  ,  for  x  h  7,  n,  £  >  0  , 

n 

»  0,  otherwise.  (A  8) 
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i/e  b 

Then  r>  =  cr  *0,11  the  true  sc&le  parameter  (e  is  also  known  as  the 
characteristic  life).  Note  that  for  the  special  case  when  0  ■  1,  the  Welbuli 
distribution  (Eq.  A7  or  Ad )  reduces  to  the  2-paraoeter  exponential  distribution 
(Eq.  A6). 

The  first  derivative  of  Eq.  A 8  with  respect  to  x  gives  the  Velbull  density 
function, 

0-1  Q 

f(x)  -  £  (2^1)  exp  [  -(JilZ_)P]  for  x  *  7,  r\,  0  >  0, 
r  * 

•  0,  otherwise.  (A  9) 

The  instantaneous  failure  rate,  Z(x)  for  the  Velbull  case  Is 
6-1  "ft 

d  M(x)/dx  ■  0(x-7)p  /ep  which  Is  a  decreasing  (increasing)  function  m  (x-7) 

If  0  <  1  (0  >  l)  and  a  constant  equal  to  l/n  if  0  -  1,  l.e.,  the  exponential 
case.  A  plot  for  f(x)  vs  (x-7)  la  shown  in  Figure  1  for  ■  1  and  for  several 
values  of  0,  the  Velbull  shape  parameter. 

The  variable  X  nay  be  standardized  by  Introducing  the  relationship 
Y  -  (X-7)/rj  .  Then  the  Velbull  cumulative  distribution  function  (c.d.f. )  and 

the  corresponding  probability  density  function  (p.d.f.)  of  Y  are  the  following: 

0(y)  ■  1  -  exp  (  -y^  ]  ,  for  y  >  0,  0  >  0, 

-  0,  otherwise.  (A  10) 

g(y)  -  0  yP_1  exp  (  -y^  ]  ,  fox  y  >  0,  0  >  0, 

-  0,  otherwise.  (A  11) 

These  are  known  as  the  standard  Velbull  distributions.  Since  Figure  1  shows 

f(x)  vs  (x-7)  for  t)  ■  1,  the  figure  also  represents  g(y)  for  the  same  set  of 
0  values. 

In  order  to  find  the  moments  of  X  It  Is  sufficient  to  calculate  the 
moments  of  Y.  The  moment  of  Y  about  zero  is 

EY*  ■  /  yk  0  y^’1  e’^dy  ■  /  zk^e"*  dz  ■  T  (£  +  l)  (A  12) 

o  o  P 
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From  this  expression  all  ocnents  of  X  may  be  found.  Hie  folloving  are  the 
nean,  standard  deviation  and  skewness  of  X,  for  b  »  l/0, 


it  -  7  +  r(b+i) 

(A  13) 

a  -  p  (  r(2b+l)  -  r*(b*-l) 

(A  lM 

a3-  [  r(3b+i)  -3  r(2b+i)  r(b+i)  +  2r3(b+i)  ]/o3 

(A  15) 

These  relationships  show  that  the  standard  deviations,  a,  and  ciy  are  free 
of  7,  and  the  mean,  u,  contains  a  linear  tem  of  7. 

In  the  case  for  which  7  is  known,  the  variable  (X-7)  will  behave  In 
accordance  with  a  Welbull  distribution  of  the  fora  given  by  Eq.  A  1 6  which 
corresponds  to  Eq.  A  8  with  the  value  of  7  equal  to  zero.  Hence,  for  this 
report  where  7  Is  assumed  to  be  known,  only  the  following  expressions  Eqs. 

A  l6  &  A  17  will  be  referred  to  as  Welbull  distributions: 

F(x)  ■  1  -  exp  f  -(x/p)^)  ,  for  x  *  0,  p,  3  >  0 

■  0,  otherwise.  (A  16) 

•  The  corresponding  p.d.f.  Is, 

f(x)  ■  (0/p)  (x/p)^”1  exp  (  -(x/p)P), 
for  x  *  0,  p,  0  >  0 

■  0,  otherwise.  (A  17) 

Figure  1  referred  to  previously,  show’s  Eq.  A  17  precisely  for  the  various 
0  values  indicated. 
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APPENDIX  B 


Graphical  Estimation  of  Welbull  Parameters 

In  this  report,  both  the  shape  parameter,.  0,  md  the  location  para¬ 
meter,  7;  are  assumed  to  be  knovr  These  assumptions  are  made  in  vlev  of 
empirical  evidence  from  experience  vith  field  and  laboratory  failure  data  that 
indicates  for  some  components  that  these  tvo  parameters  tend  to  be  constants 
over  a  vide  range  of  applications.  The  approximate  magnitude  of  these  tvo 
parameters  must  be  estimated  from  prior  data  in  order  to  use  for  acceptance 
purposes  the  sanpling  plans  reported  here.  This  appendix  outlines  tvo  slqple 
graphical  methods  for  estimating  not  only  these  tvo  parameters  but  also  the 
scale  parameter  r\  .  Since  the  tvo  methods  to  be  described  will  Involve  plotting 
the  lifelength  data  on  probability  papers,  each  of  the-  tvo  methods  also  provide 
a  graphical  test  of  fit  to  varn  against  the  use  of  Welbull  sampling  plans  if 
the  fit  is  "poor".  The  tvo  methods  I  &  II  of  estimation  are  as  follovs: 

Method  I 

This  method  is  baaed  on  the  fact  that  Eq.  Afl  can  be  vrltten  as  the 
folloving  upon  taking  the  logarithm  tvice: 

In  In  — i -  *  -  P  In  n  +  p  In  (x-7) 

1-F(x) 

*  -  In  a  +  0  in  (x-7)  (Bl) 

Hence  on  a  graph  paper  vith  In  versus  In  -  In  coordinates  the  Welbull  c.d.f. 
will  appear  as  a  straight  line.  The  folloving  sketch  depicts  such  a  straight 
line. 


In  In  — : - 

l-F(x) 
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The  Weibull  probability  paper  (Figure  4)  has  four  scales  at  its  borders,  two 

horizontal  and  two  vertical.  The  upper  and  right  scales  are  called  the  principal 

scales,  The  principle  scales  are  linear  because  they  are  for  the  plotting 

of  ln(x-7)  versus  In  In  — L__  which  are  linearly  related  according  to  Eq.  Bl. 

l-F(x)  ^ 

The  lover  and  left  scales,  the  auxilllary  scales,  are  non-linear  being  calibrated 

for  direct  plotting  of  raw  data.  The  two  additional  vertical  scales  along  the 

rlghthand  border  of  the  paper  labeled  \i/r\  and  c/r\  are  calibrated  for  0  up  to 

7.0  in  accordance  with  Eqs.  A13  and  Al4,  i.e.,  u/t>  ■  r(h+l)  ignoring  7  and 

c/n  ■  (r(2b+l)  -  r2  (b+l)]^.  These  scales  facilitate  calculation  of  |i  and  a 

as  will  be  illustrated  later.  The  additional  scales  that  appear  at  the  top  of 

the  graph  paper  are  for  snail  values  of  0  (0.02  to  1.5). 

Consider  the  following  hypothetical  lifelength  data  (in  units  of  100  hrs.) 

consisting  of  ten  ordered  observations,  x^  i  x^  i  ...  i  x^q,  for  which  we 

estimate  F(x1)  by  the  so-called  plotting  positions:  i/(lO+l),  for 
i  «  1,2,..., 10.  (See  the  last  paragraph  of  this  appendix  for  other  plotting 
conventions.)  The  observations  are  ■  27. 5*  31*  34,  38*  4l,  44*  47,  51*  57  and 
64.  The  corresponding  plotting  positions  are  ■  .09,  .18,  .27*  .36,  .46,  .55* 

.64,  .73,  .82  and  .91.  These  ten  points  plotted  on  the  Weibull  probability 
paper  give  Curve  A  (Fig.  4)  which  upon  extending  toward  the  bottom  scale 
gives  7  ■  15-  If*  however,  the  smallest  observation,  Xj  •=  27-5*  were  chosen 
as  an  estimate  of  7,  the  data  after  subtracting  27-5  from  each  observation 
will  be  *=  0,  3.5,  6.5,  10.5,  13.5*  16.5,  19-5*  23.5,  29-5  and  36-5.  With 
plotting  positions  remaining  unchanged,  the  Weibull  plot  vill  then  be  Curve  B. 

The  fact  that  these  tvo  curves  (A  and  B)  have  opposite  curvatures  indicates 
that  the  true  location  parameter  7  lies  somewhere  between  15  and  27*  5 -  Upon 

A 

several  trials,  7  «  20  (or  thereabouts)  gives  the  adjusted  observations  as 
7.5*  11*  14,  18,  21,  24,  27*  31*  37  and  44.  Using  the  same  plotting  positions 
as  before,  these  adjusted  observations  give  Curve  C  which  is  approximately 

AN  ^ 

linear.  Curve  C  gives  In  a  e  .96*  O  ■  390  (x  in  hundreds  of  heirs),  0  ■  I.85, 

A 

n  -  27-5  x  100  hours.  Some  other  life  quality  measures  may  now  be  estimated 
by  a  few  simple  calculations,  thus  (in  hours),  (Refer  to  Figure  4), 

Mean  ■  2000  +  .888(2750)  •  4442  (See  Eq.  A13) 

Std.  dev.  -  .5(2750)  -  1375  (See  Eq.  A  14) 

Reliability  function  at  2300  or  R(2300  -  2000) 

-  R(300)  -  1  -  .019)  -  .981 
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Reliable  life  at  90*  -  2000  +  780  -  2j80 
Reliable  life  at  95*  -  2000  +  530  »  2530 
Median  life  -  2000  +  2230  -  4230 
Initial  failure  rate  per  100  bra.  *  .26*. 

Since  (x-7)  instead  of  x  are  plotted,  the  initial  failure  rate  referred  to 
in  the  last  entry  is  actually  the  failure  rate  of  components  at  age  »  2000  hr 8. 
The  Curves  labeled  L(y^)  and  U(y^)  are  lower  and  upper  confidence  bands  on 

F(x)  vi-.v  c<  nfidence  coefficient  •*  80*.  For  calculations  to  determine 

21 

these  curves,  readers  are  referred  to  Kao 
Method  II 

Ibis  method  is  based  on  the  relationship  Y  »  (X-7)/r»  Introduced  in 
obtaining  Eq.  A10  from  Eq.  AS  which  is  shown  in  the  following  sketch. 


1 

P' 

0 


-0(y)  ■  1  -  a 


-/ 


-(— > 
i  1  i  ' 


o 


x 


From  this  process  of  standardisation,  we  see  that  for  0  <  p*  <1, 
p’  *  F(x)  ■  0(y).  Therefore  we  may  write, 

y  ■  0“1Cp *  J  ■  0”1(f(x)]  (B2) 

IHrt  the  relationship  Y  ■  (x-7)/rj  is  the  same  as, 

X  «  qy  ♦  7  (B  3) 

which  is  the  equation  of  a  straight  line.  Hence  on  a  graph  paper  with  the 
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Y-axis  calibrated  in  accordance  with  Eq.  B2 ,  the  Weibull  c.d.f.  wll]  again 
be  a  straight  line.  As  a  witter  of  fact,  the  above  argument  holds  for  all 
distributions,  Weibull  or  otherwise,  where  F  and  G  are  the  original  and  the 
standardised  c.d.f.  respectively.  Unfortunately,  for  the  Weibull  case,  the 
standard  c.d.f.,  G,  and  hence  also  its  Inverse  function  G~\  depends  on  0. 

(This  is  not  so  for  the  normal  distribution  or  the  exponential  distribution.) 

As  a  result,  for  each  value  of  0,  it  is  necessary  to  give  a  different  calibration 
for  y  from  Eq.  B2.  Ibis  is  precisely  what  is  done  for  this  second  kind  of 
Weibull  probability  paper  (Fig.  5)*  She  y- scale  calibration  for  values  other 
than  integer  values  of  0  may  be  interpolated  between  indicated  Integers. 

For  0- value s  of  1.0  and  3*0  (Curve  D  &  E),  the  reversal  of  curvatures  was 

A 

noted.  A  few  trials  with  other  values  resulted  in  Curve  F  with  0  ■  I.85. 

A 

For  0  *  I.85  it  was  found  that  the  plot  of  original  data  was  approximately 
linear.  Extending  Curve  F  toward. the  bottom  border  gives  7  •  2000.  The 
estimate  for  the  scale  parameter  n  >  2750  is  obtained  by  taking  the  difference 

A 

between  the  X -value  corresponding  to  Y  ■  1  and  7  which  is  the  x- value  correspond- 

A  ^ 

ing  to  Y  ■  0.  Note  that  r*  is  the  slope  of  the  line  labeled  Curve  F. 

There  are  many  plotting  conventions  for  estimating  p^  »  F(x1)  in  addition 
to  i/(n+l)  used  here.  The  following  gives  a  listing  of  some  more  widely  used 
ones: 

(1)  Sample  c.d.f.,  i/n 

(2)  Symmetrical  Sample  c.d.f.,  (i-f)/n 

(3)  Mean  of  c.d.f.,  i/(n+l) 

(4)  Mode  of  c.d.f.,  (i-l)/(n-l) 

(5)  Median  of  c.d.f.,  H"1  (I  |  i,  n-i  +  l) 

(6)  C.d.f.  of  EX1,  G(EY1)  «  F(EX4) 

(7)  Corrected  c.d.f,  ( 1-3/8 )/(n+i-) 

The  merits  of  some  of  these  plotting  conventions  are  discussed  in  Chernoff  and 

22  21 
Lleberman  for  the  normal  distribution  and  in  Kimball  J  for  the  normal  and 

24 

Type  I  extreme-value  distribution  as  well  as  in  the  text  by  Blom  which 
concerns  itself  mainly  with  the  properties  of  plotting  convention  (7)  in  the 
list  above  for  underlying  normal  distributions. 
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